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ABSTRACT 


The torsion flutter and forced response of tuned and mistuned 
cascades is examined using a standing wave approach, as opposed 
to the traditional traveling wave approach used in cascades aerodynamic 
models. The motion of the blades and the corresponding cascade aero- 
dynamic loads are expressed in terms of standing wave modes and arbitrary 
transient motion, by fitting the sinusoidal force coefficients in terms 
of ratios of polynomials in the Laplace transform variable, sometimes 
referred to as Pade approximants . Whitehead’s two dimensional, incom- 
pressible aerodynamic model is expressed in this transient form and is 
used to solve the flutter and forced response problems. Results obtained 
with the transient, standing wave analysis for the flutter and forced 
response are similar to those obtained by traveling wave analyses, but 
they yield the transient decay rate associated with vibrations of the 
blades, as opposed to the structural damping required for flutter obtained 
by the traveling method. 


The standing wave analysis presented here may prove to be more 
versatile for dealing with certain applications such as mistuned rotors, 
"localized" blade flutter, low engine order excitation, transient impulses 
on the rotor, and coupling in with forced response and dynamic shaft prob- 
lems. 
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NOMENCLATURE 


^A ? A 8? A& jAj 


A + ,A' 


AK,A t , A r 7 Aj 
A, ft 


A,B,C,t> 

a- 

CL 


^tfi'fiefitfi'fio 

-Pj. * B I A ^0 


d r c £ # 

p fr p *?^j } g H* 


nondimensional traveling wave force coefficients 
as they enter (2-40) and (2-41), 

^ A AR +iA AX^ etc * 

nondimensional aerodynamic coefficients 

associated with forward and backward 

/ + + \ 

traveling waves, * ( A AR +iA A*)» etC# 
defined by (2-10) 

matrices of first order form of equations of 
motion, defined by (3-13) for the tuned 
cascade, and by (3-33) for the mistuned 
cascade 

coefficients of cubic equation, (3-8) 
real part of eigenvalue 

elastic axis location in semichords aft of 
midchord 

coefficients of transient, standing wave 
form of airloads 

matrices of coefficients of transient, standing 
wave form of airloads, defined by (2-38) for 
the tuned cascade, and by (2-42) for the 
mistuned cascade 
semichord 


force and moment coefficients as defined 


l o 
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CCk) 
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t 



h 

h 

I 
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by Whitehead [ 2 ] 

Theodorsen Function 
matrix defined by (3-17) 
matrix defined by (3-31 ) 
imaginary part of Theodorsen Function 
vector of applied loads in first order form 
of equations of motion 
nondimensional coefficients associated with 
components of engine order excitation 
real part of Theodorsen Function 
coefficients of transient form of airloads 
matrices of coefficients of transient, standing 
wave form of airloads, defined by (2-38) for 
the tuned cascade, and by (2-42) for the 
mistuned cascade 
displacement of blade 
magnitude of blade displacement 
identity matrix 
mass moment of inertia 
average mass moment of inertia 
reference mass moment of inertia 
inertia associated with the high and low 

frequency blades of a mistuned cascade, 
respectively 
blade number 

reduced frequency, cob/U 
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¥ 



torsional stiffness of blade j 

iL 3 

Laplace transform 

Jl 

lift on section of isolated airfoil 

rrt 

aerodynamic moment acting on section of 
isolated airfoil 


aerodynamic moment acting on section of 
blade j 

D 

IY1» 

J 

disturbance moment acting on blade i 


D 

m. 

nondimensional disturbance moment acting on 
blade j 

N 

number of blades 

n 

engine order 

P 

transformation matrix of blade coordinates 
and multiblade coordinates 

P 

nondimensional eigenvalue, s <*./*>,+ T u>/ u) 0 

%cr f *lsr 

multiblade coordinates 

t 

vector of multiblade coordinates 

I 

/v 

eigenvector of multiblade coordinates 

r 

mode number, number of nodal diameters 

s' 

Laplace transform variable 

R 

radius of rotor to section of blade under 
consideration 

T 

matrix defined by (4-7) 

t 

time 

Li 

velocity at blade 

U. 

reduced velocity, U./u3 0 b 





vector of modal coordinates of first order 
form of equations of motion 
coordinates associated with augmented 

equation of standing wave form of airloads 
vector of coordinates, Y . Y 

9 cr 9 A *r- 


rotation of isolated airfoil 
rotation of blade j 

complex magnitude of torsional vibration 
modal coordinate associated with traveling 
wave having interblade phase angle 
magnitude of torsional vibration of 
blade j 

interblade phase angle, - 2ifr/N 
undamped natural frequency ratio, 
frequency ratio of vibration, m u)/u> 0 
Kronecker delta 
structural damping 
aerodynamic damping 

elastic axis location in % chord aft of 
leading edge 

position of blade j with respect to fixed, 
non-rotating axes 

position of blade j with respect to axes 
fixed to disk 


* 
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-P 

T 

S 

T 

y 

to 


XL 


inertia ratio 

stagger angle 

density of air 

nondimensional time 

phase angle of vibration of blade j 

Kussner Function 

frequency of vibration 

natural torsional frequency of vibration 

of blade j 

natural frequencies of high and low frequency 
blades of a mistuned cascade 
reference frequency, average frequency 
angular velocity of rotor 


Subscripts 


t 


i 

s 

1 

j 
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) 

matrix 

/V 


♦ 


: ) 

d/dt 

O 

) 

d/d'C 


3 


\ 


•i 

i 


blade number 


1 . Introduction 


Host work involving the aeroelastic stability and response of 
cascades makes use of cascade aerodynamic models which assume that all 
blade motions are sinusoidal. Modes of motion are represented by waves 
traveling around the circumference of the rotor. Dugundji [l] has shown 
that the traveling wave modes can be expressed in terms of standing wave 
modes, which are traditionally used in solving the fixed wing flutter 
problem. One advantage of this method is that the equations of motion for 
rotating structures and "static" structures can be easily joined. 
Eventually, problems of blade-disk coupling and shaft motions can be more 
easily handled by this method. The casting of the cascade traveling wave 
airloads into standing wave form is detailed in Chapter 2 . 

Whitehead [2] gives a good history of the early development of 
cascade aerodynamic models. His model for a two-dimensional, 
incompressible fluid with flat plate airfoils is used in this work. Smith 
[3] has extended to subsonic compressible flow, and Adamczyck and 
Goldstein [4], among others, to the supersonic cascade with subsonic axial 
flow. More recently, computational fluid dynamics has been used by 
Whitehead [5] to obtain the airforces for the two-dimensional high 
deflection cascade with subsonic axial flow. In each of these models, the 
force coefficients can be expressed as a complex number depending on the 
frequency of vibration and other flight and geometric parameters. Using 
these methods to generate airload coefficients requires a lengthy 
computation for each frequency of vibration desired. To solve the flutter 
problem the V-g method is used, whereby the structural damping required 
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for flutter is obtained. 

When control end aeroelnntlc lnternctions of fixed »ing airornft are 
hein* examined, the airloads must he in a form suitable for general 
motions of the airfoil. Edwards, Ashley and Breahuall [*] have shown that 
through the use of Laplace Transforms, the complete unsteady airloads 
acting on an isolated airfoil can be approximated by a ratio of 
polynomials, and that a transient form of the airloads can be expressed. 
The present work extends these ideas to oasosds aerodynamics, expressing 
the airloads in a transient form so that the airloads due to general 
motions of the airfoils can be determined. In saving the flutter 
P-Oblem, it is no longer necessary to know the frequency of vibration «hen 
computing coefficients. In addition, the actual transient decay of the 

“*i“" 18 0bUi " ,d - 1,1 C1 “P‘” 3. the transient form of the airlo.ds is 
applied to solving the torsion flutter problem. The isolated airfoil 

stability is compared to the tuned cascade stability, and the effect of 
mistuning is introduced. 

Kistuning, or differences in natural frequencies, stiffnesses, and 
inertias batmen blades is being studied ns a possible passive control of 
turbomaohin.ry flutter and forced response. Srinivasa. [ 7 ] ha s 

demonstrated the beneficial effect of mistuning on torsion flutter snd 
Was. and Ki.lb [8] have shoe, that mistuning often ha. a beneficial effect 
on the bending- torsion flutter speed and adverse effect on the forced 
response of the cascade. The present sort duplicate, these reeult. for 

the torsion degree of freedom ueing airforce, expressed i„ . transient, 
standing wavs form. 
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Ths effects of mistuning on the forced response are discussed in 
Chapter 4. The fact that the airloads are expressed in terras of general 
motion allows a quick solution of the response problem. If the airloads 
were expressed in terms of traveling waves, then the airloads must be 
recalculated for each frequency. To demonstrate the transient form of the 
airloads, the response of a cascade to an impulsive loading of one Made 
is calculated. 

Chapter 5 deals with the extensions of these standing wave methods 
to deal with flexible disk rotors. Also these methods are applied to 
give the cascade aerodynamic influence coefficients, that is, the effect 
of the motions of one blade on the forces produced at another blade. 
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2 . Transient Cascade Airforces 

Unsteady cascade aerodynamic theories have been under development 
since the early 1950 s. Each of these theories makes the assumption that 
the airfoil motion is sinusoidal. Edwards et. al. [6], and Vepa [9] 
give a method for obtaining airloads due to general motions of an airfoil 
using the Theodorsen function. This method was developed to study active 
control of aeroelastic structures, and is based on the approximation of 
the unsteady air loads by ratios of polynomials in the Laplace transform 
variable, sometimes refered to as Pade approximants. 


2.1 Standing Wave Representation 

The cascade aerodynamic theories were developed for harmonic 
vibrations that can be represented by traveling waves moving around the 
circumference of the disk. For a tuned rotor, each of the blades will 
have the same amplitude and interblade phase angle for a given mode of 
vibration. The torsional amplitude, <rtj, of each blade in such a cascade 
can be written in terms of a traveling wave mode as, 

_ r(<ct+j/3 r ) 

°<j = K r« (2.-1) 

- llLC » interblade phase angle 


where 
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- blade number 


N = number of blades 


= mode number, number of nodal 
diameters 


- modal amplitude 

See figure 1 • 


The airforces corresponding to such a vibration can be expressed as, 

2 . 1 , \ c(uit + 4/3) 

ori ( 2 ,- 2 ) 

Here the aerodynamic moment per unit span due to pitching motion, m^. , of 
each blade is given in terms of a traveling wave mode, where the 
nondimensional complex aerodynamic coefficient ( A^+iAj) is associated 
with the interblade phase angle /3 r . For a mistuned cascade, the blades 
no longer vibrate with the same amplitude or constant interblade phase 
angle. In order to specify the N magnitudes and N phase angles of the 
motion of each blade of an arbitrarily mistuned N-bladed rotor, the 
solution must contain 2xN constants. To obtain solutions, then, one suras 
over all nodal diameter traveling wave modes, 




m 


0 
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ORIGINAL PAGE W 
OF POOR QUALITY 



= Z ** r e 

TiO 


U-3) 


Here, there are N modal amplitudes, <* r , and N modal phase angles, /3 . 
When summed, the blade amplitudes and blade phase angles are 
reconstructed. The corresponding airforces are summed also as follows, 

tr).*ZirpU- b £ ( A R + :A i) r cK r e (Z-4) 

J r-z o 


An alternative, based on standing wave modes is discussed by Dugundji 
[l]. The use of standing wave analysis is traditional when studying the 
flutter of aircraft wings. In the present analysis, deflections are 
represented by standing wave modes of the form 

(2.-5) 

where ^. r is a generalized coordinate, and is in general a function of 
time. For example, if 0 = 0 the motion represented would be the 
umbrella" mode, where all blades move in unison. To show that the 
traveling wave representation and the standing wave representation are 
equivalent, add two traveling waves of the form (2-1), traveling in 
opposite directions (+j and -j) with equal amplitudes and take the 

real part, 
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ORIGINAL PAGE <9 
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cos (ut-+jp r ) +• 



Z 


coS 


Lwt-jf %) 


12 -b) 


Expanding and combining terms leaves 

o( = S< r coscot cosj/3 r 



which is of the form (2-5) for |Jt)=*co3<ot. The form of the standing 
wave airforces can be determined in a similar manner. Adding the 
traveling wave airforces associated with +j and -j directions gives, 



( 2.-20 


H e re, (A** +iA + ) refers to the complex aerodynamic coefficients associated 
with the traveling wave having an interblade phase angle /3 r , while the 
coefficient (A>ilC) is associated with the traveling wave having an 
interblade phase angle -j3 r , or (560-/3,.) . By examining the motion 
associated with these traveling waves, the direction of the wave can be 
determined. Interblade phase angles between /3 r - 0 degrees and /3 r - 180 
degrees are associated with forward (in the direction of rotation) 
traveling waves, while interblade phase angles between /? r - -180 degrees 
and p> r - 0 degrees are associated with backward traveling waves. 
Expanding (2-8), taking only the real part and combining terms, yields 


ORIGINAL PAGE 19 
OF POOR QUALITY 


~ ZiryUL b $[a r °( r co&ait — Aj Ccs^fi r 

+[- A r ^ c osu)t - ^^ r SvA(*itJ 

( 21 - 9 ) 

where 

A* = z ( A R, + A r ) 

^1 Z ( Ax + A x ) 

“ 2l (A r - A r ) 

^1-iCA^-Aj) (z-| 0 ; 

As has already been shown, by comparing (2-7) and (2-5) the expression for 
the torsional displacement of blade j can be written in terms of an 
arbitrary function of time and a standing wave mode as, 

*i M \ ( 2-11 ) 

The corresponding standing wave form of the moment acting on blade j can 
be written as. 
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mj » 2i tflff {[ B, \ | cr + B e £ |c(3 Cos i/ s r 

+ (2-12) 


Here the aerodynamic moment is given as a function of the displacement and 

velocity of the generalized coordinate, 9 . By assuming the solution 
. b< f 

^ cr =0 ( t cosoot, and comparing (2-9) and (2-12), one can see that the 
standing wave coefficents B,, B e , , B o , are related to the traveling 
wave coefficients A R , A t , A R , Aj, by the following relations: 

B 0 = A R 

k B, = Aj 

K---% 

kB, = A r ( 2 -i*) 


Since the traveling wave coefficients A R , A^, X^, and A t vary with 

reduced frequency k - <ob/U, the standing wave coefficients B , B^ , B ( , 

A* 

and B e can be used to fit the traveling wave data over a range of 
frequencies. Equations (2-15) can be solved for the standing wave 
coefficients. The general motion form of the standing wave aerodynamic 
moment (2-12) contains information about how the moment would vary for 
different frequencies of sinusoidal motion, as well as for any arbitrary 
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motion. The generalised coordinate <^is an arbitrary function of tine, 
igain, tne atanding wave coefficients are in general a function of the 
interblade phase angle, /3 . g.p-to-chord ratio, a/c, and stagger angle, 

f ' * S “ eXmpU ' coe TTicients fro. Whitehead's [ 2 ] incompressible, 
two-dimensional theory were fit using this si.ple approximation. The 
approximation to the Whitehead moment using (2-12) is shown by the dashed 
line in figure 2. I„ this example, the stagger angle is aero degrees and 
the forward and backward traveling wave coefficients are identical. Thus, 

V and l z are each equal to aero. The remaining standing wave 
coefficients are related to the Whitehead coefficient £ by 

a r + 1 = 2-C h * fe-14) 

»s can be seen in figure 2 this fit ( 2 - 1 5 ) might be acceptable for low 
values of reduced frequency, but is poor at higher reduced frequencies. 
Tor a more accurate fit of the data the moment acting on blade ) can be 
more generally expressed as, 



- Y 

a 


C r 


3.X 


= iq 

cr U 


cr 


(2-lfc) 
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This representation was motivated by the complete unsteady airforces 
acting on an isolated airfoil (see [lO]). If the coefficient G = 0, 
then B represents a virtual modal mass, B t a modal aerodynamic 
damping, and B 0 a modal aerodynamic stiffness. The coefficient G, 
represents the unsteady aerodynamic effect, or the amount of aerodynamic 
lag, since by solving (2-16) for Y^. one obtains 

ft *” (j — *r\ 

^ t k U-n) 

For a value of at time T" , where T is any time back to 0, the 

effect on the moment (2-15) increases as t is taken closer and closer to 

X • In the limit as t » t , the term G,Y r (t) reduces to G.| <f (t) , 

another term of the form of an aerodynamic stiffness. That is, the 
farther in the past the motion <5 occurred, the less effect there is on 

Per 

the aerodynamic moment. Hence, an aerodynamic "memory" or "lag" effect. 

In general, G, and tf, contribute to the stiffness and damping, so that 

B, and B„ alone do not represent the level of modal aerodynamic 

damping and stiffness, respectively. Note that only one augmented 

equation is used. Another possible augmented equation for Y with q is 

<r Jc 

neglected. This does not result in a significant error and it reduces the 
number of equations describing the moment acting on the blades. Thus for 
each second order equation expressing the aerodynamic moment for general 
motions of the airfoil, there is only one first order augmented equation. 
Solving (2-16) by Laplace transform and substituting into (2-15) gives 




pn 
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£> ["j] = 2-ITJ5 u x b^[ s% B, s + 


B. 


+ ?A -*- £ % 5 % bJ 


st'/ o 


*f* 6. 4* 

C> 


Wl 


(£-ie; 


where s = e'b/u, and a' is the Laplace transform variable 
to time, t. The unsteady effect is contained in the PadS 

G l 3 / 9+ 3c and 0,8/ 3+3" . 


corresponding 

approximants. 


To obtain the relations between the standing wave ooeffioients B^, 

1 3®’ '&* * and ?5 ( , and the given aerodynamic 

coefficients ij , £*,and Aj, in terms of reduced frequency, k 

solve (2-16) by substituting * cr -^oo.* t ♦ T^sin^t. Then solve for 
Wj»nd Y tar and substitute back into (2-15). Use the solution 
^-<* r cos«t in (2-15) and compare terms with (2-9) to give the new 
relations between the standing wave coefficients and the traveling wave 
coefficients as follows, 
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- k* £> + B 0 + = A r 

kB, + k<| 0 €r, s Aj 

k‘ + 9" 

1 /v /v 2 — 

-k B, + B 0 ■+ k &, = -Aj 

k*+ 3f 

kB, + k 3 0 S, = A r (2.- 19) 

There are many methods to obtain fits to the traveling wave 
coefficients using the nine standing wave coefficients B^, B ( , B 0 , G ( , 

/V AS / v " A 

g , B 2 , B J , B^ , and Gj . Since reduced frequencies at which flutter 
usually occurs are in the range k * 0 to k * 1 , that was the range chosen 
as the one over which to fit the data. To obtain the standing wave 
coefficients, the traveling wave coefficients A R , A^, and were 

fit at k * 0, k * 0.1, and k * 1, and (2-19) was used to solve for the 
standing fave coefficients. Appendix A has a more complete description of 
the fitting process. Relations (2-19) give a much better approximation to 
the Whitehead coefficients than the simpler relations (2-13). (See figure 
2.) The coefficients A^, X Jf A R , and A x are related to the forward 
and backward traveling wave coefficients by the relations (2-10). 


The solution for ■istunad rotors requires a summing over all the 


lameter .ode,. t„ addition to the ooefjft) .odee. the e.,ndin 8 

9l ” <49 " ) m,,St Sl3 ° b9 re P««»ted for a oomplete eolation. The 

standing .are demotion for ein( W , nan he obtained ei.ilarly by 

subtracting the two traveling waves «,.+ 

g waves that were added to give (2-7). 

Reducing leaves, 


— <x r siv\ wt s\t\ ^fS r 


(z-zo) 


Shioh ie juat a specie! oaee of the etanding d.flaotion pattern. 


*J U M sln dA 


(z-2i) 


^ ° f “ * P*iP of .area traveiing in oppoaite 

been ahoen equivalent to a pair of etanding n m . 

orreeponding airforce, give r.eult. of the for. <».„). (2 _ 16) . lnd 

( 2 - 19 ). 

Summarizing, to repreeant arbitrary .otion and correepondin, 
airforoee for a caeoad. naing etanding .eve .odae, the d.fleotione and 
correapondlng airforoee can be expreeaed in the folloeing fora. 


litko Sj> r+ 


(2.-ZZ) 
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°p poor 


p A QE jg 

Quality 


m i “ * VB, ? tr + B 0 j cr+ G, Y cr 

’®* U* hr" ®' u ?sr" ®» ?sr" ^ Ys3 co5 ^r 


( 2 .- 13 ) 







( 2 .- 21 ) 


(2.-25) 


where B^, B |t B #t 
obtained by fitting 
(2-19) and (2-10). 


G l ’ % ’ V 3 » 

the given sinusoidal 


, and G, , are constants 
traveling wave coefficients by 


The above procedure can be extended to both bending and torsional 
motion of the airfoils by expressing the traveling wave deflections, h f , 

and rotations, <x^ , and the corresponding traveling wave lift, , and 
moment, ra # , ae, 
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l (w*t + $Pr) 

e 



C^-aO 



c C ul ‘^ + 3/ 8 r') 

e. 


where h^ - and « f are generalized coordinates associated with the 
traveling wave mode r, and A^, A R , A,., and A^ are nondimens ional 
aerodynamic coefficients for lift due to deflection, lift due to rotation, 
etc. The coefficient A^ refers to either (A^+iA^) or (A A1L +iA At ) 
(depending on the direction of the traveling wave), as they appeared in 
(2-8). The sign conventions are defined at the leading edge as in figure 
3. The relations between the traveling wave coefficients defined above 
and Whitehead traveling wave coefficients are as follows, 


t-fcCp, 

A e = " C fM. 

A 

A g = (2 ”27) 


The transient form of the deflections and the corresponding lift and 
moment acting on blade J for a single pair of standing wave modes can be 


expressed as follow. , 


original pag c 


V’ K Ce>i iPr+ Kr^'iPr 

°<j = \ f Cos-j fa + o( sr 5 ; M ^ ^ 
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X £_ n cr p k Kr +-> U 

“ u z fc +B ,„ u T +&„ + fr^Y 


. — _ hll' — D k hsr _1> ksr — V 

Z^ u b D |AU b D 0 * b &IA 'Wsr 




&, e Y* 


r 

k iw . *> t h 




k ^2i_f • r k jl<r i o k=. r \ p* v 

\U X b + & IA U b + 4>A b ^*lA \cr 


t** ... 


U a "sr + & iB U * Sr + &oe*sr + & ieY.<*r 


S U* + ^B W *cr +■ B * + tlgY^Jsil 
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U 'Cr * 9o Y,, tv . 

u ic, - % Y,c r 

0 

0 

0 

# 

m 


The form of m*/b is similar to the form of L. By replacing all the 

b :a b y B <c » and a11 B tg by B*^ , one can obtain the expression for 

moment mj/b from the expression for lift, A*. There is one augmented 
equation (2-30) for each of the four coordinates, h , h , , and o< . 

c «* Sr <r Sf 

The equations for lift and moment can be written more compactly in matrix 
form as follows, 


U 


— 

U c r 
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Zir j> Ui b V 


{hi% 


+ B,;| + 


<? + e, Y ] 

V sx* 


H t3.Y - f 1 

^ /V ^ 


(z-3z; 


where the matrices have been defined as 


A 


- 





Y = 




Yu.. 


kc \ 


oCcv 


Y. 




Y 


*sv- ; 


( 2 .- 33 ) 


P» 


COS 


5 A 


cos 


SA 


o slnjfir J 



Ut 1 
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B, = 


B,e 1 

A/ 

z& 

B,_ 1 

- _l 

/V 

/V 

-®zc 

/V 



— 

| 

iV 

S 2 A 


| 


3 

2T> 


3 .= 


Th« matrices B, , B„ and 0, ar, aimllar to B,. 

ThU »ork involvae only tha toralon dagraa of fraadom. and only tha 
ooafflcient U uaad diractly. Tha oonatanto ra,uir.d aW. .an. 
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obtained "sing Whitehead's [ 2 ] two-dimensional, incompressible 
aerodynamics. Equations ( 2 -, 9 )) .are used to mol,, for the standing .aye 
coefficients, and their values are listed in Tables , and 2 for stagger 
angles of - 0 degrees and J - 45 degrees and gap-to-chord ratio, 
3/c-K The values listed in the table refer to an axis system located 

*t Uadi °® ed8e - Fieures 4 5 she. the fits of the Whitehead 

traveling wave coefficients at the leading edge. Tables 3 and 4 list the 

eiact values of the Whitehead coefficients versus the approximate values. 
To transfer information to any axis, located a distance t] per cent 

chord behind th. leading edge, one finds th. following relations. (See 
figure 3 .) 


COoj = (a a ) c 

- 2t) (A a ) # 

(aJvj = (A c ) o + Zrj (A A \ 

(A 0 ) 7 - ( A„) # + 2r)[(A t> ) t + Ca c ) ] + t^A.) 


(2-3+) 


.here ( )_ represents a quantity evaluated at the leading edge. The 
same transformations are used when transform* standing wave coefficients 
^0 any axis T| . By substituting B # , B^, etc. into 


the above for 
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etc. one can obtain the appropriate transferred coefficients. 


Returning to the cascade torsion problem, expressions (2-22), (2-23), 
(2-24), and (2-25) can be written more compactly in matrix form as follows 



( 2-33 ) 


[h&i +p, u i 
k ui + 3o r 


+e. r ^ 


(2.-34.) 

(? - 37 ) 


where the matrices have been defined as 


t * L tosjft. sinjA r J 





( 2 .- 36 ) 


The matrices B ( , B 0 , and G, are similar to B , and all the entries 
of all the matrices depend on the interblade phase angle. To represent 
mistuned cascade aerodynamics it is necessary to include all nodal 
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where, for r = 0, there is no contribution from the sinj/S mode, and so 
it is neglected. Since the standing waves appear in pairs and the sinj^. 
mode is a non-mode, odd numbers of blades are easier to handle. To model 
an even number of blades, an extra set of modes would be required. The 
mode number r would be N/2 , and the sinjyS r mode is always zero, since 
N is an even number. Only one mode, the cosjft. mode, representing 180 
degrees of phase separation of the motion of adjacent blades, is added. 
For simplicity, only odd numbers of blades are used in this work. The 
matrices B, , , and G, , are similar to B. . 

2.2 Application to Isolated Airfoil 

For comparison, the complete unsteady airforces acting on a typical 
section of an isolated airfoil have been expressed in terms of the general 
motion of standing wave modes, and the standing wave coefficients B t , 

B i * B c * » and . have been determined. These airforces are 
given by 


I 
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+ C(k)(*H) [-ii+«+(i-a)i^]| (2 - 


where C(lc) is the Theodorsen function (derived in terms of sinusoidal 
motion), and a is the elastic axis location in semichords aft of 
midchord. See reference [lO]. To obtain the traveling wave coefficients 
in the proper form, assume solutions of the form 

h = K a. cu>t f 

( 2 ' +s ) 


Now the airforces (2-43) and (2-44) can be represented by 
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% = IrrpLLkl A a - + (2-4f>; 

T)l/b= 2 TTj>lLb^ t "L 4 /\ D »< ^ e. (2.-47) 


where for these isolated airforces the nondimensional aerodynamic 
coefficients have been defined as follows. 


A e = W 2 ’+ t i W + CU)[j+ 

A t = ^ak - i k ( Q + i) C 

* i ( i + q *’J k*- i + (a^^)C(K)[l+Lk(z-A)] 


(Z-4&) 


If the Theodorsen function is written C(k)-F(k)+iG(k) , then we find the 
following 

A**» ik%k& 

A At * -kF 

Aba= i«k%F_(±-<OkG- 

A w » ik + 6 + Ci-0 KF . 

Acs 1 + (.«+£) k& 

Aci»-C<+i)kF 

A,«- {(. i +«*) k a +(o + i) F M'* 1 ) kG- 
A ot = - z(.^-a)k +(«i£)& kF 


* 1 * ff|tr 
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Now these coefficients can be fit using the . 

g the first two equations of (2-1 9) 

as the Theodoraen function C( k ) - » . » u . hnown Action of reduced 

~ ir - ~ * solved for fr0m (2 ., 9) „ 

shown in Appendix A. The values of fh. i 

values of the transient coefficients ere shown 

a «U 5 for . . (elastic axis at ,h. l.adin, .d,.,. Por 

' airfoil, the nodal dlaplacenent, are „u.l to the hlads 

displacements, and ( 2 - 28 ) becomes, 


hj - h (t) 
<*(t) 


(z-50) 


The corresponding transient lift and moment are obtained from ( 2 -29) as 
follows , 

f i = 2ir f U b ? B 2A^ f + « 4 -*■ B 0A t + G- lA Y h 

+ ** * + B .b« A + V + S-.bY. ] 

tyWey eft { B tc £ £ ^^ + Bot h + ^ 

+ 8 *»&«+ +<5^^ 

+ 3.)' h 8 u h 

( 2 .- 51 ) 

where the modal cress-coupling term* ... 

upllng tern, do not sppesr for the isolated 


r 


-*F**PH 
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airfoil case. The approximate aerodynamic coefficients are compared to 
the Whitehead values in figure 8, and their numerical values for several 
reduced frequencies are compared in table 6. 

Notice that the fits to the isolated airfoil aerodynamic coefficients 
are not quite as accurate as the fits to the cascade aerodynamic 
coefficients. Rock and BeBra [ll] have shown that for an isolated airfoil 
in a wind tunnel, the accuracy of the approximation depends on the ratio 
of the wall spacing to the semichord of the airfoil. They have indicated 
that to keep the same level of accuracy the order of the approximation 
must be increased as the wall-spacing-to-semichord ratio increases. This 
implies that for a given order of approximation, the aerodynamics for a 
cascade with a gap-to-chord ratio of one will be more accurately 
represented than these of an isolated airfoil, which has an effective 
gap-to-chord ratio of infinity. 




Application to Torsion Fluttei 


3.1 Isolated Airfoil 


written 


Por an isolated airfoil, the torsional equation of motion can he 


* = m+m° ( 3 -,) 

where m is the aerodynamic moment and «? is a disturbance moment acting 
on a section of the airfoil. Sign conventions at the leading edge are 
shown in figure 3. Prom the fitting of the isolated airfoil unsteady 
airforces, the moment m can be expressed in the form of (2-15) and (2-16) 


2 

m = 2f/ u * * 3 ^ +Sib xJ (3-2) 


b V 




( 3 - 3 ) 


where the standing wave coefficients shown in Table 3 have been 
transferred to the midchord are as follows , 
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= --°' 5 » W,..* >«* , 
= - 500 ^ (O. = --' s 8, 

n — • * 


3 .- 


= , /3 O 


Slno. there 1. only one de 8 ree of freedom, there ie a single multiblade 
coordinate ehich is equal to the rotation of the airfoil (i.e., 1 - «). 

By placing (5-2) into (3-1) and dividing hy ltyW 4 «£, the nondimensional 
equation of motion can be obtained 


4 +V£* =2B„S? +2CB 1P 5 

2 o ^ 


+ 2u\« +2u’'& tI> Y (3-4) 


V + U a Y = 


(3-S) 


•here the following nondimenaional quantities have been defined, 
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^ TTfV* 1 U 'fc = ' C , )=( ) t 

V- “s£ 77 JL- 

V o) t > a= « 0 b /•■*-£ 


To look at the flutter problem, n? baa been set to aero. The aubaoript 

O. denoting .a moment eoefficient due to blade angle of attack, .in be 
dropped at this time. 


To study the aeroelastic stability. 


assume a solution of the form 


— 

e 




(3-7) 


where 



Substituting ( 3 - 7 ) into (3-5), solving for r, substituting 
(3-4), and rearranging gives a cubic equation in p 


for 


Y 


in 



4 Bp* 


4Cp + JD - 


o 


/ n N 





0 
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t 

f 

i 



A = -V ~ 2.B Z 

B= Vg e i; 4-2?'1))C-2 B z 9 <> m -zb,u ( ,. 9 j 

C= Zu-SKj^+ull-Z (B.3.+ B.+ G-,) d* 


The cubic equation is solved for the roots 3 = a /«_> + i jf , and a /«?. . 

Flutter occurs when p = ±i* and a x ± 0. The roots of this equation are 
shown by the dashed line in figure 7* Only the upper half-plane is shown* 

The plot shows the roots associated with increasing values of reduced 
velocity u, , for a structural damping ratio "5=0, and frequency ratio 
^ = 1 . Using data for NASA Test Rotor 12, the value of the inertia ratio 
“V ■ 86.2, for the elastic axes of the blades at their midchords. The 
addition of an augmented state produces a pole with a high decay rate. 

This pole stays on the negative real axis as the reduced velocity is 

; 

increased, and is not plotted. • 

As can be seen in figure 7, the isolated airfoil does not flutter for 
this set of parameters. At a reduced speed of zero the free vibration 
response is obtained. Since there is no structural damping, the free 
vibration is an undamped sinusoidal vibration. As the reduced velocity, 

U. , is increased, the exponential decay rate (the magnitude of which is 
given by the real part of the eigenvalue) of the motion increases due to 
aerodynamic damping only. Finally, as the frequency (the imaginary part j 

of the eigenvalue) drops to zero, the airfoil becomes overdamped and no 


€ 
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*ibr.tory response i. possible. Hots that the root loos, is shown only up 
to - 4 , at which speed the response is still a damped vibration. 
After the root loco, touches the real axis, it moves in the direction of 
the positive real axis. The point of divergence is reached when the root 

in at the origin. Increasing the reduced velocity further brings the root 
into the unstable right half-plane. 


3.2 Tuned Cascade 


For a tuned cascade the equation of motion fo/ any blade is identical 
to the equation of motion for any other blade. Here we look at the case 
of a tuned cascade, where hy definition no structural coupling of the 
blades is allowed. This is also the case of blades on a rotor with a 
massive, rigid disk. The equation for each blade can be written 



C* . -f 

J 


2 I 

) J 



D 


3 I 


(3 -/£>) 


8y substituting for the aerodynamic moment as obtained 
chapter (2-36), .nd dividing by , one 
nondimensional equations of motion: 


in the previous 
can obtain the 
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v iV 


J rv /w ^ * 

+ £ %+TL,^ ^ V ? + nru (*5-//; 

~ ~ ~ J j 



(3'iz] 


where the nondimensional disturbance moment acting on the blade is 

I »* rrv / . 

= / 1Tj>W <«£~ 


and the matrices are defined by (2-38) n ov + 

«« oy u 3a). Next, express the blade rotatioi 

on the right hand side of (3-1 H -in +^™» * 

U ) terms of multiblade coordinates < 

mode r as, 



(2.-29) 


Premultiplying by 



yields, 



}-<0 
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i r w % + f i = 

/v ~ 

= p t p{^b 4 0 { + 2.uB, 

^2-U* & Y 1 + P T «*? (3-i*) 

+ MJ. Y = J (3-' + ) 


Rearranging yields, 

£1*4 + C 9 + K ? 


* 

+g-y =p T|w r (3-/s^ 
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W * f T Vjp - Zf T p B, 

^ ^ ^ ^ /v /w 

S ’r^w&r ? 5- 

K = P T V, *,* f -2 «» B c c-,7 

G = -25 2 p T f &, 

~ ~ ~ 

B * *3° 

r* 

This set of two second order equations (5-15), and two first order 
augmented equations (3-16) can be cast into six first order equations of 
motion for the tuned cascade as follows, 

O 

AX - B X *F 


fvis’j 
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To look at the flutter problem, let P - 0 and aaaume solutlono of the 
form 

2 ! ~ X e ^ C5-20) 

mhioh give, upon substitution into (3-18) and rearranging, the standard 
eigenvalue problem 

S I =r A B X (3-2/) 

Ip— ** 0 * ^ ' 


Tne real and imaginary part of the eigenvalues fS - a/o>„± itf , r , ploMed 
in figure 7 for the case of a nine-bladed cascade. Plotter occurs uhen 
the real part of the eigenvalue. a/ % . 1 . mere, corresponding to a 
solution of the form, 



49 


pm is 

0F POOR QUALfTY 


X = X 


— * 


This is an undamped sinusoidal motion that characterizes flutter. If the 
frequency ratio, ^ , drops to zero before the exponential damping rate, 
a/co c , the phenomenon of divergence is observed, as the damping rate 
crosses the imaginary axis. The solution in that case grows 

exponentially, 

I = X a.>o 

The system shown in figure 7 is a nine-bladed rotor with a stagger angle 
y = 45 degrees, gap-to-chord ratio of 1, zero structural damping, an 
inertia ratio ■ / \} » 86.2, and the elastic axis at midchord, rj » 0.5. The 
plot shows the eigenvalues associated with each interblade phase angle 
mode. For each interblade phase angle there are two roots, one 
corresponding to a forward traveling wave, the other corresponding to a 
backward traveling wave. The three modes that flutter are forward 
traveling waves. There will be more about the flutter modes in a 

following section. For a tuned cascade with nine blades and a stagger 
angle of 45 degrees, the flutter reduced speeds associated with an 
interblade phase angle mode are shown in figure 7. These were verified by 
comparing with tho sinusoidal solution results using Whitehead 
coefficients directly. Sec Appendix B. A comparison of the structural 
damping required for flutter, and the actual damping for a given reduced 
velocity are shown in figure 8. The difference between the sinusoidal 
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solution curve and the curve ^ =* 0 is proportional to the amount of 
structural damping required to obtain flutter. The curve obtained by 
solutions of the form J e ?t shows the decay rate associated with 
vibration in that mode for a given reduced velocity. Since there is no 
structural damping, the decay rate indicates the level of aerodynamic 
damping present. The value of reduced frequency k - u>b/u at which 

" 0 i3 the flutter reduced frequency. The two dampings, although 
essentially different in origin, are approximately equal. P or the 

important region below the flutter point the two values are hardly 
distinguishable . 


3.3 Mistuned Cascade 

MUtuning refer, to small differ.oces between the natural 
frequencies, atiffnesses, and inertias of the bladea of a oascada. Since 
it is easiest to measure blade frequencies in actual rotor, the 
undamped natural frequenciss of the blades are the quantities to be 
mistuned. For an, blade ) the natural frequency is defined as, 

W *9 = (3 -2.Z) 

•here refers to the torsional spring stiffness of the blade and I . 
is the moment of inertia of th. bl.d. about the elastic amis. » reference 


frequency is defined as the average natural frequency of the blades as 
follows, 

± V 

If one blade is to have a 5 % higher frequency than the remaining blades, 
the frequency of the blades must be given by, 



or, lolving the above to give the blade frequencies, 

|,o5 N N 

6) . . s — ' to (0 CO 

K, u v ' N+- .05 0 UvJ V , 0 S’ • 

Which yields, for a nine-bladed rotor, N - 9, with one blade at a 5 % 
higher frequency than the remaining eight, 

w v. 5 i 1 “'‘'>^4io 0 = 0,99^ w„ 

To achieve these frequencies one must mistune the blade stiffnesses and/or 
biis blade inertias. If frequency is to be mistuned by varying the inertia 
alone, then the stiffness is constant and (3-22) can be used tc solve for 


the inertias 
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J low 
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\o W 


( 3 -'.'+) 


This means that the average inertia is not equal to the reference inertia. 


L AVfr _ 


1 

M 




For a nine-bladed cascade with one blade of 5% higher frequency than the 
remaining blades, the difference between the reference inertia and the 
average inertia is slight, I o « 1.00/6 

For alternately mistuned blades, the difference between adjacent 
blades is set to, say, 5*. To find the high and low frequency for an 
average (reference) frequency of 1 , 



(3-25) 


or, solving for the frequencies, 

2 > t M 

^•jW 3 2.«S fl + .os 


ZNUe 

2.05N+.0S 


For N ■ 9 with 5$ "alternate" mistuning, 
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^ 0/373 

The corree ponding average inertia, I - I.OOI3 I o 

Tor a raistuned caac.de the equation of motion of each Made i, a, , or 

caacade (3-10). To expresa the motion of each blade of a 
miatuned caacade, the modal coordinatea are anmmed aa follow, 





This allows for N different blade amplitudes and N different blad 
P nglea. The moment acting on a miatuned caacade (3-15), (3-16), va 
developed in the previona chapter. By aubatitntlng for the blad, 
amplitude, and moment into (3-10), and rearranging a. „ae don. for th, 
tuned caacade, aimilar equation, of motion are defined aa folio.., 


+ + = p tn* (3-Z7.) 

00 m 00 mm 0 m 00 0m 

f+HY'=i ( 5 - 28 ) 


Tor a miatuned cede 1th S blad.., there are H ..oond order equation, 
of motion (3-27), and B oorre.pondin* augmented equation. (3-28). The 


' «4n A r. 
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new, N-by-N matrices have been defined as follows, 

n- p t N 3 - Jf - zp t £ 

~ z-' '' z- •**' ~ 

K=f T [v,i.^f -ZuVfO. 

& ® — Z u P P G , 

H= ag. 

All the matrices needed above are defined b.- (2-42). By making 
orthoganality of the standing wave modes, one obtains, 

N-* 2 V-/ ^ 

£ COS jfi r * Z sin jfi T = /z. 

AH 

• ^ 

3=0 


(3-2.9 


use of the 


(3 - 30 ,) 


Sr applying (5-30), th. matrix product f T f , yhlch appaara In (3-29) oan 
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be reduced to a very simple form. 


P t P = d 

~ r» ) 


D = 


N 


% 


0 


% 


P*' = D' ' P T 


o 


where the first entry is N because it represents the sum 
H-\ 

H 6t>S -1*0 * Kl 

0 J 

Casting (3-27) and (3-28) into first order form yields 

Ai - BX = F 

where the matrices of the 3N order system have been defined as, 






(3-51) 


(3-3Z.) 
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A flutter mode will now consist of a linear combination of all interblade 
phase angle modes. Setting the disturbance vector P * 0, the 6-bv-6 
standard eigenvalue problem (3-21), now becomes a 3N-by-3N eigenvalue 
problem: 

p I - A 1 B X 

Because of the form of the 3N-by-3N matrix k, only the N-by-N matrix M 
must be inverted to set up the eigenvalue problem. 

A computer program was written to set up and solve the eigenvalue 
problem for the tuned (3-21), or mistuned (3-34) cascade. See Appendix C. 
The root locus associated with the system with or.e blade of a nine-bladed 
rotor being mistuned 5% in frequency by varying the mass is shown in 
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figure 9- Each of the root loci of the mistuned system approach the root 
loci of the tuned system for sufficiently large values of reduced 
velocity, it • The relationship between reduced frequency, ^ , and 
reduced velocity is, 


r. tub w b 2L 

K 'T : u,T !! u (5-35) 

All roots on the plot with a value of u greater than about 1 are within 
the range of the original approximation to the unsteady airloads, k = 0 
to k * 1. Note the effect on the flutter speed. For the case of one 
blade mistuned 5%, the flutter reduced speed increased 1956. The flutter 
reduced frequency decreased by 1856. The magnitude of this effect 
decreases with an increasing number of blades. The one mistuned blade 
becomes a smaller part of the system and localized behavior is observed. 
To show that this is the case, a similar five-bladed cascade gained 3356 in 
flutter reduced speed, and lost 26% in flutter reduced frequency over the 
tuned case. For a realistic number of blades, one can expect that the 
flutter speed will not be increased greatly. 

The eigenvalues associated with the alternately mistuned system are 
shown in figure 10. The eigenvalues are plotted in the complex plane, for 
increasing values of reduced velocity, u . Again, for high values of 
reduced velocity the individual mistuned root loci approach the tuned root 
loci. A measure of the effectiveness of this type of mistuning is the 6156 
increase in flutter reduced velocity corresponding to a 4056 decrease in 
the flutter reduced frequency. 


The effect on the eigenvalues of mistuning one blade at a given 
reduced speed is shown in figure 11. At this value of reduced speed two 
of the tuned eigenvalues are in the unstable right half-plane. Mistuning 
just one blade shifts all of the eigenvalues. There is now one high 
frequency eigenvalue, and the remaining eigenvalues have shifted towards 
each other. This shifting together .s beneficial in that the unstable 
eigenvalues have been made more stable. The shifting has also resulted in 
less damping for the most stable eigenvalues, a problem which can have a 
adverse effect on the forced response of the system, as will be seen in 
the next chapter. See reference [8]. 

"Alternate" mistuning has a much greater effect on the eigenvalues 
(See figure 11). Now there is a high frequency group and a low frequency 
group of eigenvalues. All the eigenvalues of the mistuned system are in 
the stable left half-plane. Again, the least stable eigenvalues are made 
more stable, but the most stable eigenvalues are made less stable. The 
effect that the level of alternate mistuning has on the flutter reduced 
speed is shown in figure 12. Starting with a tuned rotor (alternate 
mistuning of 0^), the level of mistuning was increased gradually to 1 0>6. 
There is a significant effect on flutter speed, as was shown in [8]. The 
example is a nine-bladed rotor, with a stagger angle of 45 degrees, 
gap-to-chord ratio of 1.0, no structural damping, and with the elastic 
axis of the blades at the midchord. 
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3.4 Flutter Vibration Modes 


The information in the preceding sections was obtained entirely from 


the eigenvalues of tuned and mistuned cascades. To study the flutter 


vibration modes of the blades, it is necessary to examine the eigenvectors 
associated with the eigenvalue of the form f = ± itf . Solutions at the 


flutter point are of the form 


„ * £ 
i = 1 C 




where J is the complex eigenvector denoting the notion of the modal 
coordinate at the flutter point. To examine the motion of the bladea use 


(2-39) to obtain, 


2 - P f e £ * T =B ^* r 


(i - 3 ?; 


where « is the eigenvector of complex amplitudes of the blades. Motion 

of the type (3-37) can also be represented in terms of amplitude and phase 
as follows, 


= 


/5|.e‘ rarz+ ^ 


CS-38) 


where the amplitude M and phase 0 of each blade are given by 
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(- 3 - 1 . 9 ) 


To examine the actual motion of each blade, 


take the real part 


of 


(3-38) 


to obtain. 



( 3 ' 4 0) 


If adjacent blades have a different phase, then the blade vith a larger 
value of phase Is the blade that leads the motion. Figure 13 shorn, the 
amplitude and phase of each blade, for a tuned cascade, and the amplitude 
and phase for the too types of mistuning discussed in this chapter. For 
the tuned cascade all the blades have the same amplitude and the same 
interhlade phase angle. Plotted as straight lines are the phases each 
blade mould have Kith respect to blade sero for flutter in pure fi, 40 
degrees and fi, - 80 degrees interhlade phase angle traveling uave modes. 
The actual tuned cascade fluttered in a pure fi r ■ 40 degrees toward 
traveling vave mode, in the same direction as the rotation of the rotor. 
To an observer on the disk, the traveling wave speed is » /r, where w 
is the frequency of vibration. To an observer away from the disk, the 
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wave would appear to be traveling with a speed (J2.+ u>/r), greater than the 
speed of rotation, JL • 

It can be 3een in figure 13 that mistuning one blade creates a "dead 
spot" in the area of the mistuned blade. The mistuned blade and the 
blades immediately after it (those numbered higher) do not participate in 
the flutter mode as much as those farther along. The blade immediately 
before the mistuned blade vibrates with the greatest amplitude. It is 
interesting to note that the dominant interblade phase angle of the 
flutter motion is now 80 degrees, whereas the flutter interblade phase 
angle of the tuned cascade was 40 degrees. This can be seen in figure 9 
where it appears that the root locus associated primarily with the - 80 
degrees interblade phase angle mode i 3 the first to go unstable. A 
similar case was examined in which the one mistuned blade had a 5 % lower 
frequency than the remaining blades. There was no significant change in 
flutter speed between the two mistuned systems, and the mistuned blade 

created a dead spot just as the mistuning with one high frequency blade 
had* 


In the alternately mistuned cascade, the high frequency blades 
vibrate with the greatest amplitude. Pairs of blades vibrate in phase, 
with the low frequency blades lagging slightly behind the high frequency 
blades. The dominant interblade phase angle is the same as for the tuned 
cascade, namely, 40 degrees. The results are similar to those of 
reference [a]. 

The amplitudes of the complex modal eigenvectors give an indication 


of the amounts of each standing wave mode present in the flutter 
vibration. Table 8 shows the amounts of each mode present in the flutter 
mode for the tuned and two mistuned cases studied. Again, the tuned and 
alternately mistuned flutter modes vibrate predominantly with an 
interblade phase angle of 40 degrees (one nodal diameter traveling wave), 
while the interblade phase angle for the mistuning of a single blade is 
predominantly 80 degrees (two nodal diameter traveling wave). 
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4. Application to Forced Response 


4.1 Engine Order Excitation 

Consider the case of a typical section of a single blade of the 
obstructed rotor shown in figure 14. Due to the obstruction, the blade is 
subjected to a sinusoidal variation in velocity of n ■ 2 cycles per 
revolution, where n is sometimes refered to as the engine order of the 
excitation. The velocity at which the blade approaches the sinusoidal 
gusts is JLR, where R is the radius from the center of the rotor to the 
section of the blade under consideration, and SL is the speed of rotation 
of the rotor, in radians/sec. For a stagger angle -f? = 45 degrees, 
J2-R“ ® ee fig ur ® 1» 

In order to assume quasi-steady conditions, the cyclic velocity 
variations must be slow enough to allow the moment acting on the blade to 
reach most of its steady-state value. If, instead of a sinusoidal 
variation in velocity there is a step change in velocity causing a step 
change in angle of attack, the moment acting on the blade about the 
midchord can be written 

A [ 2 

where 'Y (s) is the KUssner gust function, and s-Ut/b is the distance the 
blade travels in semichords. See [lO], for example. The Kitssner function 
is shown in figure 15. The blade must travel about 20 semichords for the 
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moment acting on the airfoil to reach most of its steady-state value. If 
this unsteady effect is ignored, the error in the moment produced by a 
step change in velocity (angle of attack) will be negligible at the end of 
20 semichords of travel. 

Now, if the sinusoidal gusts are broken up into step changes in 
velocity, and the airfoil is allowed to travel 20 semichords after each 
step change, the quasi-steady moment and the unsteady moment will be 
approximately equal. Ignoring the time lag between the moment and the 
velocity will not introduce significant errors. See figure 15. For 
nearly sinusoidal velocity variations, the error between the unsteady and 
quasi-steady moments will be less than that for the step changes 
considered above. For accuracy, any blade should travel 20 semi chords for 
each half-cycle of velocity variation, 


(vdocify ) X (£•' mt p-tr Mf-Cfd t)±zo b 

JiR X J?FT - i° t (4-2; 


For the quasi- steady assumption to be valid, the ratio of the radius of 
the rotor to the semichord of a blade must be greater than a number which 
is proportional to the engine order of the disturbance. For the typical 
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section »odel uned here, the reference radius is at the three-,uarters 
span of the blade. Using data for the NASA Test Hotor 12, R/b is 
approximately 11. Using (4-2) the minimum accepteble R/b is 6.4 for a one 
engine order excitation, and 12.7 for an „ = 2 engine order. The „ork 
that folio., makes use of the ,uasi-stead, assumption, necessarily 
limiting itself to lo. engine order excitations, examples of which are 
inlet distortion or the presence of struts or supports. Whitehead [2] has 
a provision for the unsteady effect of sinusoidal upstream disturbances, 
but the, are not used in this work. To study the effect of higher engine 

order disturbances, the corresponding Whitehead coefficients would have to 
be included. 

To incorporate an engine order forcing into the equations of motion 
already developed, assume that the nondimens ional moment acting on blade 
j due to inlet velocity distortion can be expressed as, 

D* 

m j (<)= toS ' n ®i + ^sn Sin n (4-3; 


Sse[l]. Here, n represent, the engine order and §,• represents the 
position of blade j with respect to axes fixed in space, as shown in 
figure 14. The position of blade j with respect to fixed axes can be 
represented in terms of the rotation rate and an angle that gives the 
position of the blades in the rotating frame as follows, 


6 i = e j 


( 4 - 4 } 
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where = j2-rr/N, is the position of blade j with respect to an axis 
system fixed to the disk. Substituting (4-4) into (4-3), making the 
trigonometric expansion, and rearranging yields the force on blade j due 
to engine order n, 


m • = [ f en cos nJW - f S n sift nlrtl cos n&j 

+ [f sinnilti- ? sn COS njl-t] sift nSj- (4'5) 


Writing the force on each blade in matrix form yields 


m 1 « T 


' 1 


$ cn Cos r\Sli — s\v\ r\Jbi 




■fen n at -fsn Cos n£lt ) (4-0 


where 


T = 

MX! 


?££ Q 

C63 n N v 

cos n ^ J 


/v 


* 

3«V\ n Jf0 




m 


^'7) 


Equation (4-6) can be further reduced by making use of the fact that 


# 


*•?.:# flf 

of poc'fk ". ! .rrY 


r 

(^1 ^® 5 n 2.^ " f* sin n Sit 

.t. Sin nilt +f. ,co»n$l+ 


_ J(M l»JL* 


hi 


where 


f, * S £B n-if 

t> cn u *5r\ 


r+-s) 


Ho. the forcing on the blades due to an engine order disturbance can be 
written in matrix form as. 


*“■ * r 


(4-9) 


Substituting (4-9) into (3-27) give, an aquation of notion for sinusoidal 
forcing as follows. 




■ft e l '"«. T 


r+HY = 5 


(4*lo) 


where the matrices M, C. K G H j v . . „ 

•v a, ’ i» "* ™» i and i have been defined in 

(3-29) and (2-42). Since the problem is linear, solutions for different 


I® 


_ra**s**»»i 
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engine orders can be superposed. I , is only necessary to oolve for a 
single engine order response at a time. Later, constructing the response 
to any engine order forcing pattern is a matter of superposing the 
responses to the different engine orders. 

Next, examine the form of the matrix product P T . For engine 
orders n up to the maximum structural mode (here, (N-1 )/2 ), the 
trigonometric relations (3-30) can be used to obtain 
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where is the Kronecker delta, and n is the engine order chosen. 

So, in a tuned rotor, forcing in engine order n will excite only the 
n = r structural modes for harmonics r less than N/2. In a mistuned 

rotor, the cross-coupling of the structural modes generally leads to an 
excitation of all structural modes, even though a single engine order in 
being forced. In reality, engine orders higher than the number of 
structural modes used often are of interest. In these cases the relations 
(3-30) can no longer be applied. Use instead the following to determine 


Sin n & sin : 
sin r\, Gj co s r : 

cos n sin r£j = 

COS r\ & Crt r 9- = 

J 3 


2 Cn -0 ~ C05 t n + O 

•£s'nr\ Cn ^ C n f r') 

Cn-r)^ + i sinCn +0 6^- 

^ CoS (n-r)Oj + -£C©S Cn tr) ^ 


Sm PI $ • 2 0 

3*0 3 


N-I 

3L cos rr\ & ' 

j*6 7 



m * A/ ; £A/,3A/" * 
$or m 7 s A/,2/Y, 


See [l2]. Now it can be seen that the 


r a 0 structural mode is excited 


not only by the n « 0 engine order but also the n - N, 2N, 3N, etc. 
engine orders. Whenever the sum or difference of the engine order and the 
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structural .ode number equal, a multiple of the nu«ber of blades 
-111 be an excitation of that atrnctural mode by that engine order. 


there 
As an 


example, the engine orders and the structural modes of vibration they sill 
excite in a tuned rotor are listed in Tabie 9. 


the response of the blades to engine order excitations, start 
-ith the mistuned cascade equations of motion ( 5 - 32 ), where the 
disturbance moment for the engine order chosen has been defined above, 


AX -BY 

rv ^ <v 



(3 -3 z] 


The steady-state solution to 

in-Q:? 

e can be expressed as 


forcing 


that has a time variation of the form 




( 4 -/ 7 ) 


Substituting into (3-32) and solving for I yields 



( 4 - 1 + 
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Since I = 4 I P , the torsional response of the blades, 0( , can be 
obtained by extracting the modal eigenvector ^ from % and using (3-26) 


to obtain, 


2 - l °l 


Ls-Zi) 


For a tuned rotor, A, B, F, and I are as in (3-19), where m b is 
given by (4-5) • A computer program was written to set up and solve the 

forced response problem associated with tuned and raistuned cascades* See 
Appendix C. Shown in figure 16 is the frequency response of the blades of 
a tuned rotor to an n = 1 engine order excitation. It is typical of the 
response of a single degree of freedom system. 

A comparison can be made of the damping obtained by solving the 
flutter problem in transient form with the damping ratio obtained from the 

power points of the response curve. The relationships required are 


s - e fr-15) 

O w o ) N ° VvftVf' w e> 

where &/tt e ia the real part of the eigenvalue associated with interblade 
phase angle mode being forced. In this case the backward traveling wave 
is being forced. The damping ratio using the flutter data is *5 - 0.0079 
for the backward traveling wave mode. See figure 7. From the half-power 
points of the response curve, the damping ratio X » 0.0086. Remember 
that there is no structural damping, so these numbers represent levels of 
aerodynamic damping. 


72 


Also plotted is figure 16 ere the reeposses of the rotor .istuned by 
g blade zero have a 5 % higher frequency than the remaining blades. 
Boa the n - 1 engine order excitation forces all the interblade phase 

angle .odes, resulting in different levels of response for each blade. 

Per clarity only the responses of the mistnned blade and the blades 
adjacent to it are shorn in the amplitude plot. The forced response of 

adjacent blades is similar to the response of s tao degree of freedom 

system. See, for example, [l3]. The 40 degree interblade phase angle 
representing the backaard traveling aave of the forcing function is 
apparent in the phase plot. At resonances the mistuned blade, blade 0, is 
vibrating in phase aith the adjacent blade, blade 8. Instead of lagging 
the motion of blade 8 a, in the tuned case, the motion of blade 0 
leads the motion of blade 8, for frequencies betaeen their respective 

frequencies. Par aaay from resonances the tuned and mistuned 
systems shoa essentially the same response. 

The forced response of the mistuned system is aorse than the response 
of the tuned system. To explain this, examine the shifting of the 
eigenvalues of the tao systems, shoan in figure 17. The damping ratio 

aasociated aith the mode being forced decreases ahen the rotor is 
-istuned, hence the .crashing of the forced response for this particular 
case. Although the structural mode being forced is not a pure - 40 
degrees mode (because of mistuning), that mode .ill still have the largest 
response, as aill be shoan for the case of the alternately mistuned 
cascade. Par forcing frequencies aaay fro. resonance, the phase betaeen 
adjacent blades is -40 degrees, shoeing that the forcing la of the form of 
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a backward traveling wave, which is to be expected. 

The frequency response for the "alternately” mistuned rotor is shown 
in figures 18 and 19. Again, the mistuned response is worse for the 
mistuned cascade than for the tuned cascade. The eigenvalue of the 
mistuned system associated most closely with the type of forcing has been 
shifted to the right due to mistuning, resulting in lower levels of 
damping for that mode. (See figure 17.) Note that at resonances, pairs of 
high and low frequency blades vibrate in phase, while between resonances, 
the high frequency blade leads the low frequency blade in the motion. Far 

below resonance, the motion of all blades is in phase with the engine 
order forcing, 

Th. modal response of the toned and mistuned systems is shorn in 
figures 20-22. Only the amplitudes are shorn. To reoonstruct the blade 
amplitudes from the modal amplitudes, information about the relative 
phasing of the blades is required. For the tuned system, the only modal 
response is in the mode that is being forced, a, „as stated earlier. For 
the mistuned oasoade, all modes are eroited, the largest response still 
being in the mode being forced. Bote that the mode associated ,1th the 
engine order forcing has a smaller mistuned response. 


4*2 Transient Response 


The flutter and forced response of cascades can 1 
either traveling wave analysis or standing wave analj 


examined using 
s. However, the 
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important problem of th. tr.nai.nt response of cascades oan not be 
examined using traveling „av. analysis. One way to examine the transient 
behavior of oaaoadea is by expressing the motion of the blades and 
oorresponding airloads in terms of standing save modes. 


By easting the oasoade airloads into a general form. (2-40) and 

)’ the tran31ent response of the cascade can be calculated. The 
response of a tuned cascade to an impulsive motion of blade taro is shown 
in figure 23. These results vere obtained by a computer program developed 
by the Joint Computer Facility, at H.I.T. The program « (Byn a „ ic 

System Simulation) uses a fourth order Runge-Kutta integration in time. 
Blade tern undervent an impulsive motion and the subsequent "ringing down" 
of the rotor „as observed. Only the envelope of the responses is shovn in 

the lower half of figure 23. Remember that the only coupling between the 
blades is through the aerodynamics. 


With this standing wav. analysis, th. response of a cascade to 
impacts can be examined. Current transient dynamic analyse, do not allow 
for aerodynamic loads to influence th. motion of the blades. The 
technique described in this work could be applied to study the effect of 
aerodynamic loads on the transient response of rotors. 
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5. Miscellaneous Applications 


5.1 Flexible Disk Ror.ors 


The present analysis dealt with a rotor whose blades were mounted 
on a rigid disk. In many rotors, the blades would be mounted on a flex- 
ible disk so that structural coupling between the blades would be present. 
To analyze these flexible disk rotors, one generally obtains the vibration 
modes of the coupled blade-disk system by a Rayleigh-Ritz or a finite 
element method, and then one expresses the torsion angle 0^ at, s?y, the 

80% blade span section in terms of a superposition of k normal coordinates 
? (t) as, 

~ <x = # ? (s--i) 

In the above, jf is an N x k matrix relating the torsion angle o( at each 
of the N blades for each of the k vibration modes assumed in the analysis. 
For tuned rotors, the vibration modes generally occur in pairs for each 
nodal diameter r, corresponding to the Cos ^ and 5?*}^ modes, used 
earlier for the rigid disk, (2-22), (2-35). The standing wave analysis 

then expresses the rotor equations of motion in the uncoupled normal form 
as. 




where the generalized forces 

the aerodynamic moment »TJ of (2-40) and (2-41) and the 

H 

ized forces ^ in terms of the normal coordinates 
from (2-39) and (5-1) that, 



. To express 
subsequent general- 
^ , one notes 


* 



i 

i 

a 

% 

i 

| 

» 


P = § ) 


i 
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Multiplying this by the inverse of as given by (3-31) gives, 

% = f'i £ = J>''p T f } If- 4-) 

which car, be substituted into hjand then into g to obtain the general- 
ized force. 



A» 


Zn^U^b 



+ M 


U 


# 


r 

/v 


•f- 






- T, 

A- 

(f-s-J 

(5--0 


where, 

E t = PB^'V# 

/V 

E, * J t pb,jV? 

V A- ; 

Eo = f T P g.£W 

/V- A- A# A*. ^ ^ 


T, = £>Gj 

A— A-* A- /C 

3i - 3 >"p t £ 

/V ^ A ,v 

(s’- 7) 


whiie ®l ®l «•« as defined (2-42) . Equations (5-5) 

and (5-6) can be readily incorporated into the flutter and forced 
vibration of tuned and mistimed flexible bladed-disk rotors, (5-1) 


and (5-2). 


}-< 
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5.2 Aerodynamic Influence Coefficients 


For interest, it is sometimes useful to express the transient 
cascade aerodynamic influence coefficients, that is, how the motion of 
bl.ide j influences the forces on another blade i in the cascade. This 
is readily done by inverting (2-39) to give 

1 = £-o 

theu placing the above into the aerodynamic forces m of (2-40) and (2-41) 
to give. 


2 =2v,uV(£4« , E.iji T,Y} 

Y <- 3oY = X, fe. 4 (5-9J 

^ ^ A* AX 'V U ^ lO^ 

where the matrices E., , E E J J are defined in (5-7), but without 
the tor f matrices present. Allowing only motion of the j th blade 0(. and 

J 

setting all other blade motions equal to zero, gives the forces on the 

• tK# 

t blade as. 


m, = 2rT ? Ut*[( £l ). 

* £*X = 


(e,) * fe,). 0 C 3 

'ii u J v Sij i 





where Y still represents a column matrix, while T, now represents the 

* ^ 

t row of the previous square matrix £i f and X 2 the Column of the pre - 


vious square matrix The above cascade aerodynamic influence coefficients 

are valid provided the number of blades hf is large enough for cascade 
influences to die out. 
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The influence coefficients representation given by (5-11) and 


(5-12) can be written in more familiar form by solving for the aug- 


mented state variables Yfrom (5-12) and placing into (5-11) to obtain, 


- 2v ? uY{( El )..£S. ^ 


(S’- 13) 


The convolution integrals above represent a form of aerodynamic "lag" 


effect, as discussed previously, see Equation (2-17), 


Cascade aerodynamic influence coefficients were determined numeri- 


cally for the N = 3 j S/c. = 1 , ^=< 4 - 5 °, ^ = .5 case 

considered previously in most of the torsional flutter and forced response 
calculations. Table 10 gives values of the basic P matrix defined in 


(2-39) and (2-42) for a nine bladed rotor, N = 9. Table 11 shows the 


transient Whitehead coefficients used for all the torsional flutter analyses 


here, that is, type D, N = 9- S/C =1, £ = 45°, and >1= .5 cases. 

These values were interpolated from the general Table 2^ for pitching about 
the midchord (^~*5) and N =9. Finally, Table 12 gives the aerodynamic 


influence coefficients for the moment on blade i due to the motion of 


blade j = 0. For this case, (5-13) can be expressed as, 


= 2v y uV[(E 1 ). 0 ^, + 


* I (H-X. [ V (5,) " ^ it' } 


(S’- 1+) 


Because of the cyclic symmetry of the matrices, it can easily be shown 
that the due to C( 0 is the same as h^due to o(, ^ vn^due to ^ Kh 0 





■4"u 




Table 12 can be 


due to <* g , etc. Hence, 
used to assess the forces 


the coefficients given in 
on any blade due to the motion of any other 


blade. 


For comparison. 


the isolated blade (Theodorsen) fotces are also 


given in Table 12. 
a blade due to the 


These, of course, can only express the forces 
motion of the same blade. 


on 
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6. Conclusions 

The present report has given a unified standing wave approach to 
the flutter and forced response of turbine engine rotors. Both tuned 
and mistuned rotors can be readily accomodated. 

The traditional traveling wave cascade airforces have been recast 
into standing wave, arbitrary motion form, by making use of Pade approx- 
imants. Some standing wave coefficients are given for Whitehead’s 2 - 
dimensional incompressible cascade theory. 

Flutter analyses were conducted using standard constant coefficient, 
linear systems techniques. The analyses give true damping decay rates 
rather than damping margins. Typical examples for tuned and mistuned 
rotors are given. 

The forced response of the rotor to periodic engine order excitation 
of all the blades, and to the transient impulsive excitation of a single 
blade is obtained using the same aerodynamic damping as for the flutter 
cases. Typical examples for tuned and mistuned rotors are given. 

Extensions of the procedure here is indicated for the case of 
flexible disk rotors, where structural as well as aerodynamic coupling 
exists between the blades. Also, the aerodynamic influence coefficients 
for the effect of one blade's motion on the forces produced at another 
blade is given. 

The standing wave analysis methods given here can be extended to 
other flow regimes [3], [4], [5], by similar fitting of the sinusoidal 
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i. 

\ 

L- 


traveling wave coefficients. Also, they can be extended to flexible 
disk rotors as indicated previously in Chapter 5. These standing wave 
methods may prove to be more versatile for dealing with certain appli- 
cations, such as coupling flutter with forced response and dynamic 
shaft problems, transient impulses on the rotor, low engine order ex- 
citation, bearing motions, and mistuning effects in rotors. 
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Appendix A. Pitting Aerodynamic Coefficien ts 

The first step in the fitting process is to express the traveling 
wave aerodynamic forces in the form 

l -2Trj>U*bU A h + A 6 «] e ‘Cu»t + j-/9 r ) (M ) 

/b ' Zrr^U 7 1 ,[a c p +A D «J e ‘ (u,t+ iAO ( AZ ) 

where the nondimensional complex coefficients A^ , A. g , A., A , are 
functions of reduced frequency k, interblade phase angle ft , and other 
geometric and flight parameters. The relations between these coefficients 
and the tabulated data of Whitehead are given by (2-27). Traveling waves 
with interblade phase angles between 0 degrees and 180 degrees are 
refered to as forward traveling waves (in the direction of rotation), and 
those with interblade phase angles between 180 degrees and 560 degrees 
are referred to as backward traveling waves. Backward traveling waves 
have an interblade phase angle of -ft or (360 -^). m this way the 
coefficients defined by (2-10) can be generated. A^ refers to the real 
part of a traveling wave aerodynamic force or moment which has the form 


f * (A* 



t (ut+j/z ) 

« 4 e 


(A3) 


A, r.f,r, to th. real part of a traveling „,va aarodynaaic foro. or 
moment which has the form 


f =• (a: 


+ ( A j ) <k 0 e 


i( 


(A+) 





s ■ 
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Once A R , ^ A^, and A x> are formed, (2-19) can be used to solve for 

the standing wave coefficients B^, B t , B„, etc. ' Since Whitehead has 
tabulated values of the traveling wave coefficients only for k - 0, 

k ,25 ' k ~ * 5, and k s 1. the coefficients A^, A r , A R , and 
A t , were fit at k - 0, k .1, and k = 1 and the values at k = .25 and 

k = .5 were used to gauge the fit. This is a very simple way of 
determining the standing wave coefficients, and as it turns out, gives 
very good results. So, there are 12 equations; A R (k = 0, .1, 1), 
Aj(k = 0, .1, 1), A R (k - 0, .1, 1), l x (k = 0, .1, 1), and 10 unknowns, 

*’ B °’ G| * \ ’ B i» B ®» and 3c. * In order to avoid 

using two augmented equations (2-16) for each mode, and requiring two 

coefficients ^ ^ ^ was set to some intermediate value, and the 

coefficient 3, was set to zero. Some trial-and-error was involved in 

Picking the final value of ^ . There was a minor effect on the accuracy 

of the fits due to the elimination of an augmented state. The first two 

equations of (2-19) are solved for B^, B, , B p , G ( , and ^ , at k - 0, 

k = .1, and k - 1. Since the second equation of (2-19) is automatically 
satisfied for k = 0, there are 5 equations in 5 unknowns: 
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A r (0 = B 0 


|6o[A R (..)-S R( .)] = -8 2 + g -^ :()| 

- A^(o)] = -B 2 + 


10 A x Cl) • B, + 


A I (i) = B, + 


3. 6 » 
So+.oi 

9o £■ 

9. 2 + l 


Subtracting (A7) from (A6) gives 

0r f ~ /PC ^r(«0~ < ^9Ar.^°)“~A^(0 

" jT+r) 

Subtracting (A9) from (A8) gives 

q . 10 Ax («• 0 ~ At (0 

9 ‘te. “ sTrr) 


Setting (A10) equal to (All) and solving for 4 


gives 
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IP AxIO-aWi) 

100 \ (. 0 - 0 9a r O- A R CO 


(A»2) 


Adding (A6) and (A7) and solving for yields 




iooa h (.i)-ioiX k (,)+ a r (0 



(A 1 3; 


Adding (A8) and ( A9 ) and solving for B ( yields 



The standing wave coefficients B, f G, , B x , B, , can be obtained 
from equations (A5), (A12), (All), (A13), and (AH) respectively. To find 
the coefficients B, S, G, replace A R by -i £ and A r by \ 
in the above equations. The coefficients for Whitehead [2] aerodynamics 
are shown in Tables 1 and 2 for stagger angle ^ - 0 degrees and - 45 
degrees, respectively. 
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Appendix B. Flutter Points of a Tuned Cascade 

Starting with the equation of motion of any blade in a tuned cascade 








<*, +uM. 


‘3 3 


«. = rrij 


(3-0 


«here the moment n. can be written in terma of the Whitehead traveling 
wave coeffients as* 


m j = Zrj»uV[a(cJ ii ]o( e £( ' 4t ^) 


(Bl) 


Here the coefficient Cj^ is in general a function of reduced frequency, 

xnterblade phase angle, gap-to-chord ratio, stagger angle, and elastic 

axis location. See if there is a traveling wave flutter solution of the 
form 


(K ~ 


•— t C uit +j/3 r ) 

e 


Substituting (Bl) and (B2) into (3-1) yields 


(52.) 


-co ' I* +2 5H< + 1* <*£ =2ireuV [2(u<)j (S3) 

Dividing by and casting the equation into nondimensional form 

yields, 
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Where 4 , a, and '/ are defined by (3-6) and V ■ / h c , . For (B2) to be 
the solution, the real and imaginary parts of (B4) must be true 


(BS) 

Im: irij ‘ ~ [z 

(Bi) 

The real equation shows the difference between the undamped natural 
frequency ratio ^ * *0^ A*^ and the frequency ratio of vibration 
% =£0/u> 0 . Since the inertia ratio 0 is large for rotor blades in 
general, the flutter frequency is generally close to the undamped natural 
frequency, ^ ** The imaginary part of equation (B6), upon 
rearranging, gives the condition for instability as, 

^ + ^ < o 

(B 7) 

where 



(B s) 

is the aerodynamic damping. Since ^ 3^ , 

« 
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Setting the structural damping = 0 , yields the new condition for 
instability. 



f. < o 

A 


See reference [l]. The values of C M<< are tabulated by Whitehead for 
k = 0, .1, .25, .5, and 1. The approximation to the traveling rave C M<< 
can be reconstructed from the standing wave coefficients as follows. ?or 
traveling wave interblade pheje angles 0 degrees to ^? r = 180 degrees 


v- 


■\>k 2 


{b„Ic-b c , + 


k+9.* 


n 




For interblade phase angles p r 9 180 degrees to 360 degrees 


-e ! 

V W 


r ^ k + B 0 „ 




+ k Cr (D 


(Bio) 


A comparison of exact Whitehead values of *5* with reconstructed values 

A 

of ^ is shown in Table 9, for an inertia ratio 1? - 86.2. The 

example is a nine-bladed rotor with a 45 degree stagger angle, pitching 
about midchord. Only the 3 least stable interblade phase angle modes are 
shown. Since Whitehead tabulated values of coefficients at interblade 
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phase angles which are multiples of 36 degrees, the Whitehead and the 
standing wave coefficients were linearly interpolated to find values of 
the coefficients at 40, 80 and 120 degrees. 

These results can be compared with the eigenvalues obtained from the 
transient formulation. Since the tuned cascade flutter problem results in 
uncoupled equations, it is easy to look at the flutter of an individual 
mode. The uncoupled problem is third order, resulting in eigenvalues of 
the form p =£»,/#!+ [V and Since a^ is always negative and large 

compared to a ( , any response associated with a^ will be short -lived. 
The steady-state problem can be considered second order and the complex 
conjugate eigenvalues can be represented as 

^ ± l 

for example. For the tuned cascade, this means that the real 
the eigenvalue is equal to the amount of aerodynamic damping 

Qj 

OJ v J A 

Remember that there is no structural damping, and that for the tuned case 
. Now the actual aerodynamic damping obtained from the transient 
decay can be compared with the structural damping required for flutter. 
The only point whore there must be agreement is the flutter point. At the 
f- utter point, zero structural damping is required, and the real part of 


See [13], 
part of 
present , 
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> 


the eigenvalue is zero. Figure 8 shows the variation of the aerodynamic 
damping computed by both methods. 


ji 

* 






•W, 
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Appendix C. Standing Waves to Traveling Waves 

The standing wave deflections and corresponding air forces 
are given by (2-22), (2-23), (2-24) and (2-25) as. 



To reduce to traveling wave representation, one assumes harmonic 


motion of the form. 
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This implies the two standing wave coordinates anc * S«r<$ 

90° out of phase with one another. Placing the above into (c-1) gives, 

_ — liO* 

o<j = « r e s«<L]P r i- t. cs< r e 

_ iut ij'p r 


are 


= c* r e e 




= o< r e 


This gives the traveling wave deflection indicated by (2-1) . 

Placing (c-5) into the augmented state equations (c-3) 

and (c-4) and assuming harmonic motion for Y. and gives, 

J cr JSr 

, . — 'not 

(tk + 3«>)Xr e ~ ^0( r C 

— lu)t *iu?t 

(tk + <j,)X r e - t k t6f r e 


(c -0 


(c- 7 ) 


Placing (c-5) and the Xr and Ysr ^ rom (c~7) into the 


moment 


equation (c-2) gives, 


wij = 2TT5>uV|[-k l B i + Ik B, + B 0 * 5 r e lun-i P 


- [- k l B^ * i k B, + B„ + i 5, e l "<La j Pr 

+ [-k X 4 * * % + .^7| ! -] 5 r e 

(C-B) 
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Comparing the above expression with the coefficients Ad A . A© A 

K ) j K > 


defined by (2-19), one can rewrite (c-8) as. 


"Is 


= 2TT^>Ub |(A r +- i A t ) « r e 


s =? __ th)t 

+ ( Aj. “ i Ap)i(^ r 6 


- — . «_ cwt 

( A* T i A x ) to(^ 6 3 Pr 


(c-9 


+ - 


. — , = . ^ Cll)t > 

(-A £ + t A^) ot r £ 


o ..Mfr c jPr *1 , 7 . 

= 2irp Ui> j A R e -»-uA t e + A R e. -FLA^e /< 


Z7TJ> Ub | (A r + A r ) + l (Aj^ A r JjJof r C 




Then, noting the definitions of A R , Aj > A R ^ Aj in 
(2-10), the above reduces to. 


, = 2TT?Ub (a; + lAj) 5 r e 


m 


+ N » i(ut+jP r ) 


(c- 


(0 


which is the same aerodynamic force indicated by (2-2). 


Similarly, if one assumes harmonic motion of the form. 
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. . _ Lut 

%cr(*) = «, e 

%*(*) = -I <x r e 


and goes through the same reduction, one will obtain 











tf 
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»{Zb 



? 7 y , ) - , } | ; TT~r~T 


- h “ »|Z.b«< 

*rj = * 


=m + >j2.bj 


Figure 3. Force Notation 
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(Arfl+ A cg ) + l(AtI+ ^l)~ 



1 


i 


i 

Exact WV*>4eVied 
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N>^-! 

XsoW'Ve^ Aw^o'*^ 


A^oVi«H«iiev\ 

tto N ^j 
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Reduced Frequency k 


\ -X = o 

^CR ^CX 


Figure 4c. Approximate Whitehead Force Coefficients 
Type A , ^ - 0 degrees, a/o ■ 1 , *] “ 0 


I 

♦ 


X Exad vJW,+cV\eaA 































Ck 


original page m 
OF POOR QUALITY 

(a dr * A pr ) + i (a di + t BT )«2CM, 

“! I I T I 1 1 1 

X WV»-VeV\cftA 

- Ise>\qt e<J A\r£&\ \ 

■ * Aftx-fcKiMft+io* 







Reduced Frequency k 


Figure 5g. Approximate Whitehead Force Coefficients 
Type Ajj* 45 degrees, s/c - 1 , n- 
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Figure 6a. Approximate Theodorsen F v rce Coefficients 
Type A a , 
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Figure 8. Comparison of the Damping Required fo 
Flutter with the Transient Decay Rate 


* 




Real Part ^ a ' u 0 

Figure 9. Flutter Root Locus - Blade "0" Mistuned 
5 % in frequency. N ■ 9, ^ * 45 degrees 

s/c ■ 1, Tl" 0.5, “0 »86.2, 0 
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Figure 16. Frequency Response of Tuned Cascade and 
of Cascade with Blade "0" Mistuned 5 % 
to a One Engine Order Excitation 
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figure 22. Frequency Response of Nodal Amplitudes 
at an "Alternately" Nistuned Cascade 
to a One Engine Order Excitation 
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fir pP® \ 


-0.366 -0.305 0.000 0.000 0.000 0.000 0.000 0.000 0.000 

0.366 0.860 0.305 0.000 0.000 0.000 0.000 0.000 0.000 

0.366 0.110 0.000 0.000 0.000 0.000 0.000 0.000 0.000 

-0.432 -0.788 -0.110 0.000 0.000 0.000 0.000 0.000 0.000 

-0.435 -0.410 0.000 -0.055 0.140 0.000 0.000 0.000 0.000 

0.435 1.068 0.800 -0.303 0.140 0.000 0.000 0.000 0.000 

0.435 0.184 0.000 0.025 0.140 0.000 0.000 0.000 0.000 

-0.498 -0.964 -O.36O 0.138 0.140 0.000 0.000 0.000 , 0.000 

-0.496 -0.510 0.000 -0.116 0.235 0.000 0.000 0.000 0.000 

0.496 1.258 1.005 -0.319 0.235 0.000 0.000 0.000 0.000 

0.496 0.256 0.000 0.059 0.235 0.000 0.000 0.000 0.000 

-0.560 -1.128 -0.506 0.162 0.235 0.000 0.000 0.000 0.000 

-0.545 -0.590 0.000 -0.175 0.300 0.000 0.000 0.000 0.000 

0.545 1.410 1.173 -0.327 0.300 0.000 0.000 0.000 0.000 

0.545 0.316 0.000 0.094 0.300 0.000 0.000 0.000 0.000 

-0.610-1.262-0.628 0.174 0.300 0.000 0.000 0.000 0.000 


-0.575 -0.643 0.000 -0.214 0.335 0.000 0.000 

0.575 1.506 1.283 -0.328 0.335 0.000 0.000 

0.575 0.354 0.000 0.118 0.335 0.000 0.000 

-0.638-1.348-0.706 0.182 0.335 0.000 0.000 


0.000 0.000 

0.000 0.000 

0.000 0.000 

0.000 0.000 


-0.586 -0.661 0.000 -0.231 0.350 0.000 0.000 0.000 0.000 

0.586 1.540 1.321 -0.329 0.350 0.000 0.000 0.000 0.000 

0.536 0.366 0.000 0.129 0.350 0.000 0.000 0.000 0.000 

-0.648-1.378-0.734 0.186 0.350 0.000 0.000 0.000 0.000 


Tabid 1. Transient Whitehead Coefficients 
" 0 degrees, s/c - 1 , r* -0 
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Type 


B, 

B 0 

o. 

3. 

8. 

1 

** 

B , 

B o 



0.000 
0.427 
0.000 
- 1 .0081 - 0.134 


0.000 

0.000 

0.000 

0.000 


0.000 

0.000 

0.000 

0.000 


0 0.000 

0 0.000 

0 0.000 

0.000 0.000 


00 

0.000 

00 

0.000 

00 

0.000 

00 

0.000 


- 1.044 


.000 - 0.038 
.660 - 0.156 
.000 0.021 

.322 0.090 


0.000 

0.000 

0.000 

0.000 


.493 - 0.493 
.493 1.236 
.493 0.242 
.556 - 1.112 


.000 - 0.110 0.295 
.867 - 0.216 0.295 
.000 0.062 0.295 
0.452 0.118 0.295 


- 0.014 

- 0.018 

0.048 

- 0.034 


.095 
.176 
- 0.048 
- 0.068 


- 0.517 

0.517 

0.517 

- 0.580 


0.000 - 0.188 
1.041 - 0.234 
0.000 0.105 
- 0.562 0.128 


0.375 

0.375 

0.375 

0.375 


0.000 - 0.014 0.000 
0.000 - 0.015 - 0.222 
0.000 0.045 0.000 
0.000 - 0.030 0.054 


H 


.536 - 0.571 0.000 - 0.247 0.420 0.000 - 0.007 0.000 0.052 

.536 1.368 1 . 158 - 0.234 0.420 0 . 000 - 0 . 008 - 0.116 0.040 

.536 0.312 0.000 0.136 0.420 0.000 0.025 0.000 - 0.025 

.596 - 1.230 - 0.636 0.128 0.420 0.000 - 0.018 0.026 - 0.014 


0.542 

- 0.602 


.000 - 0.268 
.198 - 0.234 
.000 0.148 
.662 0.128 


0.000 0.000 

0.000 0.000 

0.000 0.000 

0.000 0.000 


0.000 0.000 

0.000 0.000 

0.000 0.000 

0.000 0.000 


Table 2. Transient Whitehead Coefficients 
■ 45 degrees, s/c * 1, - 0 


te g -ir j a 
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A- 


O.OOO-fiO.OOO 
0.004-i0.030 
0.023-10.076 
0.091 -iO. 152 
0.366-10.301 


/3 r * 


approx 


A fc - approx 


0.000+10.000 O.OOO+iO.OOO O.OOO+iO.OOO 
0.004-10.031 -0.014-10.067 -0.014-10.067 
0.023-10.076 -0.01 4-iO. 126 -0.015-10.126 
0.092-i0.153 0.059-10.220 0.058-10.219 
0.366-10.305 0.381-10.417 0.381-10.418 


A = 72 


A^- approx 


/3 r - /08 


exact A^- approx 


.0 O.OOO+iO.OOO O.OOO+iO.OOO 
.1 -0.01 3-iO. 093 -0.013-10.093 

.25 -0.031 -iO. 186 -0.031-10.185 
.5 0.029-i0.300 0.029-10.300 

1.0 0.385-10.536 0.386-10.536 


O.OOO+iO.OOO 

-0.012-10.111 

-0.038-10.233 

0.008-i0.372 

0.385-10.638 


O.OOO+iO.OOO 
0.01 2-iO. 1 1 2 
0.038-i0.234 
0.008-10.372 
0.364-i0.638 


A * 


A = ;8o 


A^- approx 


approx 


.0 

O.OOO+iO.OOO 

.1 

-0.01 2-iO. 123 

.25 

-0.041 -iO. 264 

.5 

-0.004-i0.421 

1.0 

0.383-10.707 


O.OOO+iO.OOO 
0.01 2-iO. 1 27 
0.041 -iO. 275 
0.009-10.439 
0.380-10.733 


Table 3a. Exaot Whitehead Coefficients versus 
Approximate Values, Type A ^ 

^ ■ 0 degrees, s/c • 1 , ^ - o 
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Ag- exact 


approx 




A_- e 


approx 


1 .005+10.0001 1.005+i0.000 
0.951 +i0. 010 0.951+10.011 


.25 0.803+10.155 0.805+10.155 
•5 0.619+10.507 0.620+10.506 
1.0 0.207+11.187 0.207+11.187 


14 + 


,Ol 


Ag- exact Ag- approx 


0.305+10.000 
0.301 +iO. 086 0.301+10.086 
0.282+10.216 0.282+10.215 
0.213+10.431 0.214+10.430 


0.801 +iO. 000 0.800+i0.000 

0.691-10.038 0.693- i 0.037 


1.0 |-0. 063+10. 860|-0. 061+10.860 


0.540+10.139 

0.409+10.458 


0.542+10.138 

0.410+i0.455 


0.068+11 .028 0.068+il .026 


fir * 


Ag- approx 


1 . 173+iO.OOO 
1 .134+10.042 
1.004+i0.192 
0.796+10.561 
0.329+il .320 


1.173+iO.OOO 
1.135+10.043 
1 .005+i0. 1 92 
0.796+10.561 
0.328+il .320 


1 

Ag- approx 

1.238+10.000 
1.250+10.060 
1 . 1 29+10. 21 9 
0.912+10.601 
0.412+11.408 

1.283+10.000 
1 .250+10.061 
1.130+10.219 
0.913+10.601 
0.413+11.407 j 


A " t&o 


Ag- approx 


1.290+i0.066 1.290+10.067 


Table 3b* Exact Whitehead Coefficients versus 
Approximate Values, Type A ft 
■ 0 degrees, s/c » 1 , ^ ■ 0 
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I 

jtf 

I 


% 




A- • o°__ 

exact A c - 


/<3 f - 


approx A fc - exact 


O.OOO+iO.OOO 
-0.004+10.01 i 
-O.i)??+i0.027 
-0.092+10.054 
-0.307+10.108 


O.OOO+iO.OOO 
-0.004+10.011 
-0.023+10.028| 
-0.092+10.055 | 
-0.366+10. no 


O.OOO+iO.OOO 

0.004+10.030 

-0.008+i0.056 

-0.086+10.099 

-0.413+10.188 


Pr « 72. 


A.- approx 

0.000+10.000 

0.004+10.030 

-0.008+10.057 

-0.086+10.098 

-0.410+10.187 


A c - exact 

•0 O.OOO+iO.OOO 
•1 0.004+10.047 

•25 0.000+10.093 
•5 -0.076+10.151 
1*0 -0.442+10.270 


A » iog>* 


Aj,- approx 

0.000 *- 10. 000 
0.004+10.047 
0.000+10.093 
[-0.076+10.151 
•0.440+10.269 


Pr ~ 1 44“ C 


Ag- exact 


•0 O.OOO+iO.OOO 
•1 0.004+10.068 

•25 0.006+10.145 
•5 -0.062+i0.232 
1 -0 -0.469+10 . 590 


A c - approx 

O.OOO+iO.OOO 

0.004+10.068 

0.006+10.145 

-0.062+10.232 

[-0.469+10.390 


S c - exact 

O.OOO+iO.OOO 

0.004+10.060 

0.004+10.125 

•0.068+10.199 

•0.460+10.541 


A t - approx 

O.OOO+iO.OOO 

0.004+10.060 

0.004+10.125 

-0.067+10.199 

-0.459+10.^42 


o r = 1 &£>' 


At- exact 

O.OOO+iO.OOO 

0.004+10.071 

0.007+10.152 

0.061+10.244 

•0 •472+10.407 


A«- approx 

O.OOO+iO.OOO 

0.004+10.071 

0.007+10.155 

-0.060+10.244 

1-0.471+10.406 


Table 3c. Bxaot Whitehead Coefficients versus 
Approx iaate Values, Type A c 
" 0 degrees, s/c • 1 , yj . 
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4 

\ 


A « 3.4. 


Aj- exact A, • approx A,- exact 1^- 




/Sr 

* 


V 

exact 

-0. 

.106 

-iO. 

079 

-0. 

.083 

-iO. 

197 

-0. 

,002 

-iO. 

394 

0. 

,321 

-iO. 

787 


-O.IIO+iO.OOO 


fir * Tz 


approx 


.0 -0.506+i0.000 -0.506+i0.000 
.1 -0.476-i0.055 -0.476-i0.054 
.25 -0.385-i0.202 -0.385-i0.201 
.5 -0.234-i0.503 -0.233-i0.502 
.0 0.206-il .093 0.208-i1.092 


Pr * 144-° 


Ajj- exact A^- approx 


.0 

-0.707+i0.000 

-0.706+i0.000 

.1 

-0.686-i0.035 

-0.685-i0.085 

.25 

-0.602-i0.250 

-0.601 -iO. 250 

.5 

-0.422-i0.590 

-0.421-i0.590 

1 .0 

0.096-il .293 

0.096-il . 293 


-0.361 +iO. 000 -0.360+i0.000 
-0.309-i0.031 -0.308-i0.031 
-0.225-i0.184 -0.224-i0.182 
-0.109-i0.448 -0.108-i0.446 
0.272-i0.946 0.273-i0.945 


/3 r - 


\j- exact Aj- approx 


-0.627+i0.000 -0.628+i0.000 
-0.603-i0.073 -0.605-i0.074 
-0.517-i0.229 -0.519-i0.230 
-0.346-i0.554 -0.348-i0.554 
0.141-11.214 0.142-11.214 


Pr * l^>0 


approx 


.735+ 
.714- 
.631- 
. 449- 


Table 3d. Exact Whitehead Coefficients versus 
Approximate Values, Type A^ 

^ ■ 0 degrees, s/c - 1, •'J ■ 0 
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0 , =0 


A- o c 


I A - exact 5 a - approx A a - exact A a - approx 


O.OOO+iO.OOO O.OOO+iO.OOO 
0.004-i0,043 0.005-i0.043 


.25 0.028-i0. 107 
•5 0.1 13-i0.213 


0.028-i0. 107 


0.1 13-i0.213 0.1 14-i0.214 
0.455-10.423 0.455-i0.427 


Pr *3£ 

k A a - exact A a - approx 

.0 O.OOO+iO.OOO O.OOO+iO.OOO 

.1 -0.005-i0.060 -0.005-i0.06l 

.25 0.005-i0. 1 28 0.004-i0.129 

.5 0.085-i0.235 0.084-i0.234 

.0 0.433-i0.453 0.434-10.451 


O.OOO+iO.OOO O.OOO+iO.OOO 
O.OOO+iO.OOO O.OOO+iO.OOO 
O.OOO+iO.OOO O.OOO+iO.OOO 
O.OOO+iO.OOO O.OOO+iO.OOO 
O.OOO+iO.OOO | O.OOO+iO.OOO 

fir l U’ 

A a - exact A a - approx 


O.OOO+iO.OOO 

0.028-i0.016 


O.OOO+iO.OOO 
0.028-i0.01 6 


0.030-i0.046 0.029-i0.045 
0.017-10.061 0.01 6-iO. 060 


0.434-10.451 I 0.002-i0.065 I 0 . 001 -iO . 065 


A = iz 


Pr^ll 


exact A a - approx A a - exact X A - approx 

O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO 
-0.006-i0.083 -0.006-i0.083 0.028-i0.010 0.027-i0.010 


.25 -0.01 5-iO. 177 -0.01 5-iO. 178 
•5 0.042-i0.294 0.042-i0.295 
1.0 0.392-i0.522 0.392-10.523 


0.045-i0.040 0.043-i0.040 
0.036-i0.072 0.035-i0.070 


0.392-i0.523 | 0.01 2-iO. 086 0.012-i0.087 


Table 4a. Exact Whitehead Coefficients versus 
Approximate Values, Type A a 
^ ■ 45 degrees, s/c * 1 , H ■ 0 
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ORIGINAL PAGE JQ 
OF POOR QUALITY 



fir ~ 

108* 

/3 r = I0& 

k 

exact 

A*- approx 

exact 

A^- approx 

.0 . 
.1 
• 25 
.5 
1.0 

O.OOO+iO.OOO 
-0.007-i0. 101 
-0.026-i0.221 
0.009-i0.359 
0.353-i0.600 

O.OOO+iO.OOO 
-0. 007-10.101 
-0.026-i0.222 
0.009-10.361 
0.352-i0.603 

O.OOO+iO.OOO 
0.021 -iO. 006 
0.038-i0.027 
0.037-i0.057 
0.01 5-i0.075 

O.OOO+iO.OOO 
0.021 -iO. 006 
0.038-i0.027 
0.036-i0.057 
0.01 5-i0.078 


Pr ~ 

1-+° 



k 

k 9 - exact 

X B - approx 

A^- exact 

A r - approx 

.0 
.1 
.25 
• 5 
1.0 

O.OOO+iO.OOO 
-0.008-i0. 1 1 3 
-0.031 -iO. 251 
-0.01 1-iO. 406 
0.327-i0.658 

O.OOO+iO.OOO 
-0.008-i0.1 13 
-0.031-i0.251 
-0.011-i0.407 
0.326-i0.659 

O.OOO+iO.OOO 
0.01 1-iO. 003 
0.021 -iO. 01 4 
0.022-i0.031 
0.010-i0.043 

O.OOO+iO.OOO 

0 . 01 1-iO . 003 
0.021 -iO. 01 4 
0.022-i0.030 
0.01 2-i0.044 


Pr " 1 BO° 

(5 r - i 6 

0° 

k 

A ft - exact 

A r - approx 

A e - exact 

approx 

.0 
.1 
• 25 
.5 
1.0 

O.OOO+iO.OOO 
-0.008-i0. 1 17 
-0.033-i0.26l 
-0.01 7-iO. 423 
0.318-i0.680 

O.OOO+iO.OOO 
-0.008-i0.117 
-0.033-i0.26l 
-0.01 7-iO. 424 
0.317-i0.680 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 


Table 4b. Exact Whitehead Coefficients versus 
Approximate Values, Type A 8 
■ 45 degrees, s/c - 1 , h - 0 
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O 

ll 

pr =0° 

k 

A 8 - exact 

Ag- approx 

Aj- exact 

Ag- approx 

.0 
.1 
.25 
• 5 
1.0 

0.423+i0.112 

0.399+i0.281 

0.3H+10.561 

-0.027+11.119 

0.427+10.000 
0. 422+iO. 112 
0.399+i0.280 
0.313+10.560 
-0.028+11.119 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 


fir- 

/3 r - 


k 

Aj- exact 

Aj- approx 

m 

A^- exact 

m 

Ag- approx 

.0 

.1 

.25 

.5 

1.0 

0.660+i0.000 
0.613+10.047 
0.520+i0.216 
0.400+i0.532 
0.038+il .133 

0.660+10.000 
0.617+10.049 
0.526+i0.216 
0.403+10.531 
0.038+il. 133 

0.000+i0.370 

0.112+i0.309 

0.127+i0.193 

0.082+i0.127 

0.034+i0.093 

0.000+10.370 
0.117+10.302 
0.126+i0.185 
0.079+i0. 1 25 
0.034+i0. 103 


^ o 

A- » 72 

Pr — 7 2.° 

k 

A^- exact 

A m ~ approx 

Ag- exact 

A^- approx 

.0 

.1 

.25 

.5 

1.0 

0.867+i0.000 

0.839+i0.057 

0.744+10.202 

0.582+i0.524 

0.176+11.178 

0.867+10.000 

0.840+10.058 

0.746+10.202 

0.584+10.523 

0.175+11.177 

0.000+10.308 

0.047+10.293 

0.078+10.242 

0.068+10.183 

0.030+i0.135 

0.000+i0.308 

0.052+i0.290 

0.082+10.234 

0.068+10.177 

0.030+10.146 


Table 4c. Exact Whitehead Coefficients versus 
Approximate Values, Type A^ 

* 45 degrees, s/c - 1, ■ 0 



ORIGINAL FAGE W 
OF POOR QUALITY 


lOl 


/3 r - 1 

0 

0$ 

Aj- exact 

approx 

1 .041 +i0.000 
1 .021 +i0.073 
0.937+i0.220 
0.767+i0.544 
0.319+il -237 

1 .041 +i0.000 
1 .020+i0.073 
0.937+i0.221 
0.762+i0.545 
0.319+il. 237 

A * MV 

A,- exact 

Ap- approx 

1.157+iO.OOO 

1.140+i0.084 

1.063+i0.239 

0.887+i0.568 

0.422+11.284 

1.158+iO.OOO 
1.140+i0.084 
1.063+10.239 
0.887+10.569 
0.423+il .284 

11 

CP 

Ay exact 

A g - approx 

1.198+iO.OOO 
1.181+10.088 
1 .107+i0.247 
0.931+i0.578 
0.460+il .302 

1.198+iO.OOO 
1 .181 +i0.088 
1.106+i0.246 
0.929+10.578 
0.459+11.301 


fir - lo« 


A approx 


exact 


O.OOO+iO.222 
0.020+i0.217 
0.038+10.196 
0.038+i0. 162 
0.01 7 + i0.124 


O.OOO+iO. 222 
0.022+i0. 21 6 
0.040+i0. 1 93 
0.038+i0. 161 
0.01 6+iO. 139 


/3 r * \++° 


A^- approx 


O.OOO+iO. 116 O.OOO+iO. 1 1 6 
0.007+i0.114 0.008+i0.114 
0.01 4+iO. 107 0.01 6+iO. 106 
0.01 6+iO. 092 0.01 6+iO. 098 
0.006+i0.073 0.006+i0.082 


A- - I 


O.OOO+iO. 000 O.OOO+iO. 000 
O.OOO+iO. 000 O.OOO+iO. 000 


O.OOO+iO. 000 
O.OOO+iO. 000 
O.OOO+iO. 000 


O.OOO+iO. 000 
O.OOO+iO. 000 
O.OOO+iO. 000 


Table 4d* Exact Whitehead Coefficients versus 
Approximate Values, Type Ag 
T ■ 45 degrees, s/c * 1, *) ■ 0 


m 




































143 ORIGINAL PAGE IS 

OF POOR QUALITY A 



fir**' 

/3r-0 

• 


k 

A,- exact 

A fr - approx 

A c - exact 

T - approx 


.0 
.1 
.25 
• 5 
1.0 

O.OOO+iO.OOO 
-0\005+i0.018 
-0.029+i0.046 
-0.115 + i0.092 
-0.459+i0. 182 

O.OOO+iO.OOO 

-0.005+i0.018 

-0.028+i0.046 

-0.114+i0.091 

-0.455+10.182 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

O.OOO+iO.OOO 

i 

( 


/3 r * 



k 

_ 

Ag- exact 

A t - approx 

* 

A e - exact 

A fc - approx 

i 

.0 

.1 

• 25 

• 5 
1.0 

O.OOO+iO.OOO 
0.000+i0.029 
-0.01 6+iO. 060 
-0. 100+iO. 108 
-0.451+10.206 

O.OOO+iO.OOO 

0.000+i0.029 

-0.015+10.060 

-0.099+i0.108 

-0.451+i0.206 

O.OOO+iO.OOO 

-0.009+i0.007 

-0.006+i0.020 

0.007+10.028 

0.026+10.036 

O.OOO+iO.OOO 
-0.009+i0.007 
-0.004+i0.020 
0.01 2+iO. 027 
0.037+10.029 



/3 r * 7Z° 

/3 r « 72T 

i 

i 

k 

A c - exact 

A t - approx 

A fc - exact 

A^- approx 

l 

.0 

.1 

.25 

.5 

1.0 

O.OOO+iO.OOO 
0.001 +iO. 043 
-0.005+i0.091 
-0.078+10.148 
-0.437+10.259 

O.OOO+iO.OOO 
0.001 +iO. 043 
-0.005+i0.091 
-0.077+i0. 148 
-0.436+i0.259 

O.OOO+iO.OOO 

-0.008+i0.004 

-0.009+i0.017 

0.004+i0.033 

0.033+10.047 

O.OOO+iO.OOO 
-0.010+i0.005 
-0.01 2+i0.020 
0.003+i0.036 
0.035+i0.044 

i 

i 

i 


i 


Table 4e. Exact Whitehead Coefficients versus 
Approximate Values, Type A e 

■ 45 degrees, s/c ■ 1, - 0 


i 





POOR 


qUALITt 


A - I OB 


A * \ 08 


exact 


approx 


A t - exac 


A c - approx 


.0 O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO 

.1 0.002+i0.054 0.002+10.054 -0.005+i0.002 -0.007+i0.003 

.25 0.000+i0.119 0.000+10.1 19 -0.007+i0. 01 1 -0.010+i0.014 

.5 -0>062+i0. 1 91 -0.062+i0.191 -0.002+i0.026 -0.001 +iO. 029 

1.0 -0.426+i0.317 -0.425+i0.317 0.028+i0.041 0.030+10.039 


A - 14- 4*° 


A*-? exact 


approx 


A = / 4 -V 


A t - exact A t - approx 


.0 O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO 

.1 0.002+i0.062 0.002+i0.062 -0.002+i0.001 -0.003+i0.001 

•25 0.002+i0.138 0.002+i0. 138 -0.004+i0.006 -0.005+i0.007 

•5 -0.054+i0.222 -0.054+i0.223 O.OOI+iO.OH 0.000+i0.015 

1.0 -0.419+10.359 -0.420+10.361 0.016+10.023 0.016+10.021 


A - ifio 


A t - e 


A = /&o* 


4 t - approx A c - exact A t - approx 


.0 O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO O.OOO+iO.OOO 
.1 0.002+i0.064 0.002+i0.065 O.OOO+iO.OOO O.OOO+iO.OOO 
.25 0.003+i0.144 0.003+10.144 O.OOO+iO.OOO O.OOO+iO.OOO 
.5 -0.051+10.233 -0.051+10.234 O.OOO+iO.OOO O.OOO+iO.OOO 
1.0 -0.41 6+iO. 375 -0.41 6+iO. 376 O.OOO+iO.OOO O.OOO+iO.OOO 


Table 4f* Exact Whitehead Coefficients versus 
Approximate Values, Type A t 
T - 45 degrees, s/c ■ 1, ^ « 0 
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ORIGINAL PAGE is 
0F p 0OR QUALITY 



5p- exact | A > - 


0.1 79— iO. 1 01 


0*17q + tn’m? A*^ 10,000 0 * 000+ i0-000 
0.179-i0.101 O.OOO+iO.OOO O.OOO+iO.OOO 




n 1K1 4 A OM A Ar-1 1 v.VWT 1UlUW | U.UUU+iO.OOO 

0 *am in*t?AR ~A*Acl~ i ^* 2 '’ 2 ®*°°0 + i0«000 O.OOO+iO.OOO 
: -g - • 1?*504 O.OOO+iO.OOO o.ooo+io.ooS 

— 11 ,0 ° 9 0*140-11 .008 _0.000+iQ.00Q O.nnn^n aaa 


A « 36* 


~ exact A,- approx 

Mil filpHis 

0.299-H.029| O^lKOgg -0.045-10.01 6 l-n.oZL° 



Ap- exact 


A ~ 1Z 




A* 7Z° 


exact 


-^433-iS'™ 'o'«r, l o'^ 0.000-10.089 0.000-10.088 
• 21 5-il .080 0.21 3— il .080 -0.054-10.024 -0.052-i0.*025 


Table 4g. Exact Whitehead Coefficients versus 
Approximate Values, Type Ap 
T " 45 degrees, s/c ■ 1, ^ m 0 



























oKWiiuu- 

OF POOR QUALITY 


fir ' 108 


fir * 108 


A*” exact A,- approx A,- exact A,- approx 


.0 -0.561 +iO. 000 -0.562+10.000 

.1' -0.547-i0.086 -0.548-i0.086 
.25 -0. 486-10.236 -0.486-10.256 
•5 -0.355-10. 528 -0.335-10.529 
1.0 0.15 0-11.158 0.130-11.138 


r* I 4*4-° 


0.000-i0.055 0.000-i0.054 
0.01 1 -i0.053 -0.012-10.052 
0.023-i0.045 -0.024-i0.043 
0.033-i0.033 -0.032-i0.031 
0.044-i0.022 -0.042-10. 022 


A»- exact A,- approx Ap- exact tp- approx 

.0 -0.635+i0.000 -0.636+10.000 0.000-i0.026 0.000-i0.026 
.1 -0.622-i0.094 -0.623-i0.094 -0.005-i0.025 -0.005-i0.025 
.25 -0.565-10. 251 -0.565-i0.251 -0.010-i0.023 -0.01 1-iO. 022 
.5 -0.41 1-10.551 -0.412-10.552 -0.017-10.018 -0.016-10.018 
.0 0.070-il .183 O.069-H.184 -0.024-i0.012 -0.023-i0.014 


fir * l&o' 


fir =160 


Ip- exact Ap- approx A # - exact Ap- approx 


.0 -0.661 +i0. 000 -0.662+i0.000 
.1 -0.649-10.097 -0.650-10.097 
.25 -0.592-10.256 -0.593-10.256 
.5 -0.438-10.559 -0.439-10.560 
.0 0.048-il .200 0.048-11.199 


0.000+10.000 

O.OOO+iO.OOO 


O.OOO+iO.OOO 

O.OOO+iO.OOO 


0.593-10.256 O.OOO+iO.OOO O.OOO+iO.OOO 
0.439-10.560 O.OOO+iO.OOO O.OOO+iO.OOO 
0.048-11.199 O.OOO+iO.OOO O.OOO+iO.OOO 


Table 4h* Exact Whitehead Coefficients versus 
Approximate Values, Type Aj, 

T • 45 degrees, s/c ■ 1, Q ■ 0 

































147 


ORIGINAL PAGE 19 
OF POOR QUALITY 


type 



B o 

G , 

3» 

A 

-0.470 

-0.526 

0.000 

-0.072 

0.130 

B 

0.565 

1.229 

1.000 

-0.244 

0.130 

C 

0.485 

0.263 

0.000 

0.036 

0.130 

I) 

-0.595 

-1.114 

-0.500 

0.122 

0.130 


Moments taken about the leading edge, b * 0.0, or a ■ -1 


Alternate coefficients can also ^otained using the known theoretical 
two dimensional coefficients together with a rough approximation to 
the Theodorsen function as C(p) • (.55 p + .15)/(p + .15). This 
would result in the alternate coefficients. 


type 

B 2 

B 1 

B o 

G l 

9* 

A 

-0.500 

-0.550 

0.000 

-0.068 

0.150 

B 

0.500 

1.325 

1.000 

-0.349 

0.150 

C 

0.500 

0.275 

0.000 

0.034 

0.150 

0 

-0.563 

-1.163 

-0.500 

0.174 

0.150 


Table 5. Transient Theodorsen Coefficients ^ - 0.0 




POOR 


QUALITY 


approx 


approx 


O.OOO+iO.OOO O.OOO+iO.OOO I.OOO+iO.OOO 1 .000*10.000 
-0.01 2-iO. 083 -0.022-10. 087 0.853+10.003 0.904+10.005 


.25 -0.016-10.168 -0.027-10.161 
.5 0.050-i0.299 -0.050-10. 28C 


0.737+i0.185 0.773+10.207 


0.050-i0.299 -0.050-10.280 0.586+10.548 0.630+i0.555 
0.400-10.539 0.399-i0.535 0.190+11.208 0.195+11.198 


approx 


approx 


O.OOO+iO.OOO O.OOO+iO.OOO -0.500+10.000 -0.500+10.000 
0.004+10.042 0.008+i0.044 -0.423-10.051 -0.449-10.052 
-0.006+10.084 -0.002+i0.080 -0.351-10.213 -0.367-10.229 
-0.087+i0. 150 -0.088+10.140 -0.215-i0.524 -0.237-10.527 
-0.450+10.270 -0.450+10.270 0.218-11.104 0.215-11.098 


- I ir f> U- t> l A c -g + P 


Table 6. Exact Theodor sen Coefficients versus 
Approximate Values, » a o.o 
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ORIGINAL PAGE w 
OF POOR QUALITY 


N = 9 

!?* r n 1,111 «««• th. structural 
d 5 <*»i/3 r U follow., 


°* 9, 18, 27,, 

1.8.10.17.19.26.28.. 

2.7.11.16.20.25.29.. 

3.6.12.15.21 .24.30.. 

4.5.13.14.22.23.31.. 


Talle 8. Structural Nodes Responding to Engine c 
Excitation - Tuned Rotor, H ■ 9 
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ORIGINAL PAGE 19 
OF POOR QUALITY 


*40 degrees 

(3 r - 80 degrees 

■ 120 degrees 

( 

exact 

approx 

exact 

approx 

exact 

approx 

' 1 

-0.1682 

-0.1722 

-0.1561 

-0.1588 

-0.0886 

-0.0889 

-0.0078 

-0.0083 

-0.0084 

-0.0082 

0.0008 

0.0008 

0.0011 

0.0009 

0.0013 

0.0013 

0.0034 

0.0032 

0.0011 

0.0010 

0.0012 

0.0011 

0.0015 

0.0015 


Table 9* Exact Aerodynamic Damping versus Approximate Values 
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OF POOR QUALITY 


1 

1 

0 

.1 

0 

1 

0 

1 

1 

.7660 

.6428 

.1736 

.9848 

-.5000 

.8660 

-.9397 . 

1 

.1736 

.9848 

-.9397 

.3420 

-.5000 

-.8660 

.7660 -. 

1 

-.5000 

.8660 

-.5000 

-.8660 

1 

0 

-.5000 . 

1 

-.9397 

.3420 

.7660 

-.6428 

-.5000 

.8660 

.1736 -. 

1 

-.9397 

-.3420 

.7660 

.6428 

-.5000 

-.%60 

.1736 . 

1 

-.5000 

-.8660 

-.5000 

.8660 

1 

0 

-.5000 -. 

1 

.1736 

-.9848 

-.9397 

-.3420 

-.5000 

.8660 

.7661 . 

1 

.7660 

-.6428 

.1736 

-.9848 

-.5000 

-.8660 

-.9397 -. 


See aquations (2-39) and (2-42) 


Table 10. Values of Basic P Matrix for N m 9 
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ORIGINAL PAGE 19 
OF POOR QUALITY 




B, 

Bo 

G l 

Q 

** 

B t 

/Nr 

B, 

A* 

B e 

A* 

G . 

0° 

-.069 

-.134 

.243 

0 

0 

0 

0 

0 

0 

40° 

' -.063 

-.128 

.346 

-.090 

.188 

0 

-.015 

-.240 

.188 

80* 

-.063 

-.127 

.429 

-.155 

.313 

0 

-.017 

-.209 

.144 

120 ° 

-.062 

-.124 

.493 

-.198 

.390 

0 

-.012 

-.142 

.087 

0 

160 

-.060 

-.120 

.528 

-.221 

.427 

0 

-.005 

-.050 

.030 

Isolated 

airfoil 

-.015 

-.148 

.500 

-.158 

.130 

0 

0 

0 

0 


Values interpolated from general Table 2 for pitching about midchord 
<n- .5) and for N ■ 9 


Table 11. Transient Whitehead Coefficients for Torsional Flutter 

.0 


Type D, N = 9, ^ = 45°, s/c - 1, 7} - .5 


Analysis 


























